arXiv:1504.05476vl [cs.DS] 21 Apr 2015 


Optimal parameterized algorithms for planar facility location 
problems using Voronoi diagrams* 
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Abstract 

We study a general family of facility location problems defined on planar graphs and on 
the 2-dimensional plane. In these problems, a subset of k objects has to be selected, satisfying 
certain packing (disjointness) and covering constraints. Our main result is showing that, for each 
of these problems, the time brute force algorithm of selecting k objects can be improved 

to time. The algorithm is based on an idea that was introduced recently in the design of 

geometric QPTASs, but was not yet used for exact algorithms and for planar graphs. We focus 
on the Voronoi diagram of a hypothetical solution of k objects, guess a balanced separator cycle 
of this Voronoi diagram to obtain a set that separates the solution in a balanced way, and then 
recurse on the resulting subproblems. 

The following list is an exemplary selection of concrete consequences of our main result. We 
can solve each of the following problems in time vP^'^\ where n is the total size of the input: 

• d-SCATTERED Set: find k vertices in an edge-weighted planar graph that pairwise are at 
distance at least d from each other {d is part of the input). 

• d-DOMiNATiNG Set (or (A:, d)-CENTER): find k vertices in an edge-weighted planar graph 
such that every vertex of the graph is at distance at most d from at least one selected 
vertex (d is part of the input). 

• Given a set V of connected vertex sets in a planar graph G, find k disjoint vertex sets in D. 

• Given a set V of disks in the plane (of possibly different radii), find k disjoint disks in T>. 

• Given a set V of simple polygons in the plane, find k disjoint polygons in V. 

• Given a set V of disks in the plane (of possibly different radii) and a set P of points, find k 
disks in T) that together cover the maximum number of points in V. 

• Given a set V of axis-parallel squares in the plane (of possibly different sizes) and a set P 
of points, find k squares in P that together cover the maximum number of points in P. 

It is known from previous work that, assuming the Exponential Time Hypothesis (ETH), there 
is no f(k)n°^'^^ time algorithm for any computable function / for any of these problems. 
Furthermore, we give evidence that packing problems have time algorithms for a much 

more general class of objects than covering problems have. For example, we show that, assuming 
ETH, the problem where a set P of axis-parallel rectangles and a set P of points are given, and 
the task is to select k rectangles that together cover the entire point set, does not admit an 
f(k)n°^^^ time algorithm for any computable function /. 
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1 Introduction 


Parameterized problems often become easier when restricted to planar graphs: usually significantly 
better running times can be achieved and sometimes even problems that are W[l]-hard on general 
graphs can become fixed-parameter tractable on planar graphs. In most cases, the improved running 
time involves a square root dependence on the parameter: for example, it is often of the form 
20 (Vk) Qj. ^0{Vk) ^ appearance of the square root can be usually traced back to the fact that 

a planar graph with n vertices has treewidth 0{^/n). Indeed, the theory of bidimensionality gives a 
quick explanation why problems such as Independent Set, Longest Path, Feedback Vertex 
Set, Dominating Set, or even distance-r versions of Independent Set and Dominating Set 
(for fixed r) have algorithms with running time [TTl [1211131 HU 113113 Il313 [HI SZj- 

In all these problems, there is a relation between the size of the largest grid minor and the size of the 
optimum solution, which allows us to bound the treewidth of the graph in terms of the parameter 
of the problem. More recently, subexponential parameterized algorithms have been explored also 
for problems where there is no such straightforward parameter-treewidth bound: for example, for 
Partial Vertex Cover and Dominating Set [22], ^-Internal Out-Branching and A;-Leaf 
Out-Branching IIZ!, Multiway Cut |33|, Subset TSP |3l|, Strongly Connected Steiner 
Subgraph |3, Steiner Tree [13S3I- For some of these problems, it is easy to see that they 
are fixed-parameter tractable on planar graphs, and the challenge is to make the dependence on k 
subexponential, e.g., to obtain ■ time (or perhaps • rp^^^ time) algorithms. 

Others are W[l]-hard on planar graphs, and then the challenge is to improve the known rp^^'> time 
algorithm to rP^'^^ time. For all these problems, there are matching lower bounds showing that, 
assuming the Exponential Time Hypothesis (ETH) of Impagliazzo, Paturi, and Zane |3r)[ I35j . there 
are no 2°^^^^ ■ rP^^^ time algorithms (for FPT problems) or time algorithms (for W[l]-hard 

problems). 

A similar “square root phenomenon” has been observed in the case of geometric problems: it is 
usual to see a square root in the exponent of the running time of algorithms for NP-hard problems 
defined in the 2-dimensional Euclidean plane. For example, TSP and Steiner Tree on n points 
can be solved in time jjGj- Most relevant to our paper is the fact that Independent Set 

for unit disks (given a set of n unit disks, select k of them that are pairwise disjoint) and the discrete 
fe-center problem (given a set of n points and a set of n unit disks, select k disks whose union covers 
every point) can be solved in time rp^'^'> by geometric separation theorems and shifting arguments 
I3S10E313QI, improving on the trivial rp^^^ time brute force algorithm. However, all of these 
algorithms are crucially based on a notion of area and rely on the property that all the disks have 
the same size (at least approximately). Therefore, it seems unlikely that these techniques can be 
generalized to the case when the disks can have very different radii or to planar-graph versions of 
the problem, where the notion of area is meaningless. Using similar techniques, one can obtain 
approximation schemes for these and related geometric problems, again with the limitation that the 
objects need to have (roughly) the same area. 

Very recently, a new and powerful technique emerged from a line of quasi-polynomial time 
approximation schemes (QPTAS) for geometric problems (I] [2] ESI E] • As described explicitly by 
Har-Peled [28|, the main idea is to reason about the Voronoi diagram of the k objects in the solution. 
In particular, we are trying to guess a separator consisting of 0{pk) segments that corresponds 
to a balanced separator of the Voronoi diagram. In this paper, we show how this basic idea and 
its extensions can be implemented to obtain rP^'^^ time exact algorithms for a wide family of 
geometric packing and covering problems in a uniform way. In fact, we show that the algorithms 
can be made to work in the much more general context of planar graph problems. 
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Algorithmic results. We study a general family of facility location problems for planar graphs, 
where a set of k objects has to be selected, subject to certain independence and covering constraints. 
Two archetypal problems from this family are (1) selecting k vertices of an edge-weighted planar 
graph that are at distance at least d from each other (^-Scattered Set) and (2) selecting k 
vertices of an edge-weighted planar graph such that every vertex of the graph is at distance at 
most d from some selected vertex (d-DOMiNATiNG Set); for both problems, d is a real value being 
part of the input. We show that, under very general conditions, the trivial time brute force 

algorithm can be improved to time for problems in this family. Our result is not just a 

simple consequence of bidimensionality and bounding the treewidth of the input graph. Instead, we 
focus on the Voronoi diagram of a hypothetical solution, which can be considered as a planar graph 
with 0{k) vertices. It is known that such a planar graph has a balanced separator cycle of length 
0{pk), which can be translated into a separator that breaks the instance in way suitable for using 
recursion on the resulting subproblems: each subproblem contains at most a constant fraction of 
the k objects of the solution. Of course, we do not know the Voronoi diagram of the solution and 
its balanced separator cycle. However, we argue that only rP^^) separator cycles can be potential 
candidates. Thus, by guessing one of these cycles, we define and solve rP^'^'i subproblems. The 
reason why this scheme yields an rP^'^^ time algorithm is the fact that recurrence relations of the 
form f{k) = rP^'^^ f{k/2) resolves to f{k) = rP^'^^ 

In Section [3T| we define a general facility location problem called Disjoint Network Cover¬ 


age, which contains numerous concrete problems of interest as special cases. We defer to Section 3.1 
the formal definition of the problem and the exact statement of the running time we can achieve. In 
this introduction, we discuss specific algorithmic results following from the general result. 

Informally, the input of DISJOINT Network Coverage consists of an edge-weighted planar 
graph G, a set V of objects (which are connected sets of vertices in G) and a set C of clients (which 
are vertices of G). The task is to select a set of exactly k pairwise-disjoin10 objects that maximizes 
the total number (or total prize) of the covered clients. We need to define what we mean by saying 
that an object covers a client: the input contains a radius for each object in T) and a sensitivity for 
each client in C, and the client is considered to be covered by an object if the sum of the radius and 
the sensitivity is at least the distance between the object and the client. In the special case when 
both the radius and the sensitivity are 0, this is equivalent to saying that the client is inside the 
object; when the radius is r and the sensitivity is 0, then this is equivalent to saying that the client 
is at distance at most r from the object. The objects and the clients may be equipped with weights 
and we may want to maximize/minimize the weight of the selected objects or the weight of covered 
clients. 

The first special case of the problem is when there are no clients at all: then the task is to select 
k objects that are pairwise disjoint. Our algorithm solves this problem in complete generality: the 
only condition is that each object is a connected vertex set (i.e. it induces a connected subgraph of 
G). 

Theorem 1.1 (packing connected sets). Let G be a planar graph, T> be a family of connected 
vertex sets of G, and k be an integer. In time we can find a set of k pairwise 

disjoint objects in V, if such a set exists. 

We can also solve the weighted version, where we want to select k members of T> maximizing the 
total weight. As a special case, if we define the open ball to be the set of vertices at distance 
less d from v and let V = {By d /2 I ^ ^ then Theorem |l.l| gives us an time algorithm 

^More preceisely, if the objects have different radii, then instead of requiring them to be pairwise disjoint, we 


require a technical condition called “normality,” which we define in Section 3.1 


3 





for d-SCATTERED Set, which asks for k vertices that are at distance at least d from each other (with 
d being part of the input). For unweighted graphs and fixed values of d, this problem is actually 
fixed-parameter tractable and can be solved in time using a simple application 

of bidimensionality HZ!. 

If each object in P is a single vertex and r(-) assigns a radius to each object (potentially different 
radii for different objects), then we get a natural covering problem. Thus, the following theorem is 
also a corollary of our general result. 


Theorem 1.2 (covering vertices with centers of different radii). Let G be a planar graph, 
let D,C C V{G) he two subset of vertices, letr: D ^ Z+ be a function, and k be an integer. In 
time ■ we can find a set S D of k vertices that maximizes the number of vertices 

covered in G, where a vertex u € G is covered by v G S if the distance of u and v is at most r{v). 


li D = G = V{G), r{v) = d for every v G V{G), and we are looking for a solution fully covering 
C, then we obtain as a special case d-DOMiNATiNG Set (also called (/c, d)-CENTER), that is, the 
problem of f indi ng a set S of vertices such that every other vertex is at distance at most d from 
gives an time algorithm for this problem (with d being part of the input). 


1.2 


S. Theorem 

Again, for fixed d, bidimensionality theory gives a time algorithm m- 

Theorem 1.2 can be interpreted as covering the vertices in C by very specific objects: we want 
to maximize the number of vertices of G in the union of k objects, where each object is a ball of 
radius r[v) around a center v. If we require that the selected objects of the solution are pairwise 
disjoint, then we can generalize this problem to arbitrary objects. 


Theorem 1.3 (covering vertices with independent objects). Let G he a planar graph, let T> 
be a set of connected vertex sets in G, let G CV(G) be a set of vertices, and let k be an integer. In 
time ■ we can find a set S of at most k pairwise disjoint objects in T> that maximizes 

the number of vertices of G in the union of the vertex sets in S. 


Our algorithmic results are also applicable to geometric problems. By simple reductions (see 
Section 3.2), packing and covering geometric problems can be reduced to problems on planar graphs 
(although one has to handle rounding and precision issues carefully). In particular, given a set of 
disks (of possibly different radii), the problem of selecting k pairwise disjoint disks can be reduced 
to the problem of selecting disjoint connected vertex sets in a planar graph, which can be solved 
using Theorem 0 Alternatively, the algorithmic techniques behind our main results (Voronoi 
diagrams, balanced separators, recursion) can be expressed directly in the geometric setting, yielding 
a somewhat simpler algorithm that does not have to handle some of the degeneracies appearing 
in the more general setting of planar graphs. In Section we discuss some of these geometric 
algorithms in a self-contained way, which may help the reader to understand the main ideas of the 
more technical planar graph algorithm. 


Theorem 1.4 (packing disks). Given a set T> of disks (of possibly different radii) in the plane, 
in time we can find a set of k pairwise disjoint disks, if such a set exists. 


This is a strong generalization of the results of Alber and Fiala [5|, which gives an time 

algorithm only if the ratio of the radii of the smallest and largest disks can bounded by a constant 
(in particular, if all the disks are unit disks). 

As Theorem o works for arbitrary connected sets of vertices, we can prove the analog of 
Theorem 1.4 for most reasonable sets of connected geometric objects. We do not want dwell on 
exactly what kind of geometric objects we can handle (e.g., whether the objects can have holes, how 
the boundaries are described etc.), hence we state the result only for simple (that is, non-self-crossing) 
polygons. 
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Theorem 1.5 (packing simple polygons). Given a set T> of simple polygons in the plane, in 
time • we can find a set of k polygons in V with pairwise disjoint closed interiors, if 

such a set exists. Here n is the total number of vertices of the polygons in V. 


Geometric covering problems can be also reduced to planar problems. The problem of covering 
the maximum number of points by selecting k disks from a given set T) of disks can be reduced to a 
problem on planar graphs and then Theorem 1.2 can be invoked. This reduction relies on the fact 
that covering by a disk of radius r can be expressed as being at distance at most r from the center 
of the disk, which is precisely what Theorem |1.2| is about. 


Theorem 1.6 (covering points with disks). Given a set C of points and a set T> of disks (of 
possibly different radii) in the plane, in time we can find a set of k disks in V 

maximizing the total number of points they cover in C. 


The problem of covering points with axis-parallel squares (of different sizes) can be handled 
similarly. Observe that an axis-parallel square with side length s covers a point p if and only if p is 
at distance at most s/2 from the center of the square in the l^o metric. This allows us to reduce 


the geometric problem to the planar problem solved by Theorem 1.2 


Theorem 1.7 (covering points with squares). Given a set C of points and a set T) of axis- 
parallel squares (of possibly different size) in the plane, in time we can find a set 

of k squares in T) maximizing the total number of points they cover in C. 


Hardness results. There are several lower bounds suggesting that our main algorithmic result 
is, in many aspects, best possible—both in terms of the form of the running time and the generality 
of the problem being solved. The Disjoint Network Coverage problem we define in Section [sT 
gives a very general family of problems, including many artificial problems. Therefore, it is not very 
enlightening to show that the running time we obtain for DISJOINT Network Coverage cannot 
be improved. What we really want to know is whether our algorithm gives the best possible running 


time in concrete special cases of interest, such as those in Theorems 1.1-1.7 


There have been investigations of the parameterized complexity of various geometric packing and 
covering problems, giving tight ETH-based lower bounds in many cases [3611371 155] . Many of these 
reductions can be described conveniently using the Grid Tiling problem as the source of reductions. 
Reductions using Grid Tiling involve a quadratic blow up in the parameter and therefore they 
give lower bounds with a square root in the exponent. For example, assuming ETH, the problem 
of finding k disjoint objects from a set of unit disks, or a set of axis-parallel unit squares, or a set 
of unit segments (of arbitrary directions) cannot be solved in time f{k)n°^^^ for any computable 
function / [HTIEH]. This shows the optimality of Theorems 1.4 and |1.5[ In Dominating Set 


problem for unit disks, the task is to select k of the disks such that every disk is either selected 
or intersected by a selected disk; assuming ETH, there is no f(k)n°^^^ time algorithm for this 
problem for any computable function / m- Observe that if 2? is a set of unit disks (that is, disks 
of radius 1) and C set of centers of these disks, then a subset of S' C 2? is a dominating set if and 
only if replacing every disk in S with a disk of radius 2 covers every point in C. Therefore, covering 
points with disks (of the same radius) is more general than Dominating Set for unit disks, hence 
the optimality of Theorem 1 1. 6| follows from the lower bounds for Dominating Set for unit disks. 


In a similar way, the optimality of Theorem 1.7 follows from the lower bounds on Dominating Set 
for axis-parallel unit squares m- 

As we shall see in Section 3.2, there are reductions from the geometric problems to the planar 
problems. Besides the algorithmic consequences, these reductions allow us to transfer lower bounds 
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for geometric problems to the corresp ondin g planar problems. In particular, it follows that, assuming 
ETH, the running time in Theorem 1.1 cannot be improved to (even if \'D\ = 


With a direct reduction from Grid Tiling, one can also show that there is no time 

algorithm for d-SCATTERED Set and d-DOMiNATiNG Set on planar graphs, with d being part of 
the input (these reductions will appear elsewhere). 

Comparing packing results Theorems |1.1| and |1.5| with covering results Theorems 1.2[ 1.6, and 
1.7, one can observe that our algorithm solves packing problems in much wider generality than 
covering problems. It seems that we can handle arbitrary objects in packing problems, while it is 
essential for covering problems that each object is a “ball,” that is, it is defined as a set of points 


that are at most at a certain distance from a center. (Theorem 1.3 seems to be an exception, as it 
is a covering problem with arbitrary objects, but notice that we require independent objects in the 
solution, so it is actually a packing problem as well.) We present a set of hardness results suggesting 
that this apparent difference between packing and covering problems is not a shortcoming of our 
algorithm, but it is inherent to the problem: there are very natural geometric covering problems 
where the square root phenomenon does not occur. 

Our strongest lower bound is not based not on ETH, but on the variant called Strong Exponential 
Time Hypothesis (SETH), which can be informally stated as n-variable CNE-SAT not having 
algorithms with running time (2 — e)” for any e > 0 (cf. [35]). Using a result of Patra§cu and 
Williams |11| and a simple reduction from Dominating Set, we show that if the task is to cover 
all the vertices of a planar graph G by selecting k sets from a collection D of connected vertex sets, 
then is unlikely that one can do significantly better than trying all \T)\^ possible sets of objects. 

Theorem 1.8 (covering vertices with connected sets, lower bound). Let G he a planar graph 
and letV he a set of connected vertex sets ofG. Assuming SETH, there is no f{k)-{\'D\ + \V{G)\)^~'^ 
time algorithm for any computable function f and any e > 0 that decides if there are k sets in V 
whose union covers |U(G)|. 

A similar reduction gives a lower bound for covering points with convex polygons. 

Theorem 1.9 (covering points with convex polygons, lower bound). LetV he a set of convex 
polygons and let V he a set of points in the plane. Assuming SETH, there is no f{k) • {\V\ + 
time algorithm for any computable function f and e > 0 that decides if there are k polygons in V 
that together cover V. 


The convex polygons appearing in the hardness proof of Theorem 1.9 are relatively “fat” (i.e. 


the area of each polygon is at most a constant factor smaller than the smallest enclosing disk), and 
they have an unbounded number of vertices. Therefore, it may still be possible that the square 
root phenomenon occurs for simpler polygons and we can have time algorithms. We show 

that this is not the case: we give two lower bounds for axis-parallel rectangles. The first bound is 
for “thin” rectangles (of only two types), while the second bound is for rectangles that are “almost 
squares.” 

Theorem 1.10 (covering points with thin rectangles, lower bound). Consider the problem 
of covering a set V of points by selecting k axis-parallel rectangles from a set V. Assuming ETH, 
there is no algorithm for this problem with running time f{k) ■ {\V\ -|- for any computable 

function f, even if each rectangle in V is of size 1 x k or k x 1. 


Theorem 1.11 (covering points with almost squares, lower bound). Consider the problem 
of covering a set V of points by selecting k axis-parallel rectangles from a set V. Assuming ETH, 
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for every eo > 0, there is no algorithm for this problem with running time f{k) ■ {\V\ + 

for any computable function f, even if each rectangle in T) has both width and height in the range 

[1 ~ eO) 1 + eo]- 


Theorem |1.11 shows that even a minor deviation from the setting of Theorem |1.7| makes it 
unlikely that algorithms exist. Therefore, it seems that for covering problems the existence 

of the square root phenomenon depends not on the objects being simple, or fat, or almost the same 
size, but really on the fact that the objects are defined as balls in a metric. 

Our techniques. The standard technique of bidimensionality does not seem to be applicable to 
our problems: it is not clear for any of the problems how the existence of an 0{^/k) xO{y/k) grid minor 
helps in solving the problem, hence we cannot assume that the input graph has treewidth 0{^/k). 
In more recent subexponential parameterized algorithms, we can observe a different algorithmic 
pattern: instead of trying to bound the treewidth of the input graph, we dehne a “skeleton graph” 
describing the structure of the solution and use planarity of the skeleton to bound its treewidth. 
Then the fact that this skeleton graph has treewidth 0{Vk) (and, in particular, has balanced 
separators of size 0{Vk)) can be used to solve the problem. The right choice of the skeleton graph 
can be highly nonobvious: for example, for Multiway Cut [33], the skeleton graph is the union of 
the dual of the solution with a minimum Steiner tree, while for Subset TSP |34j . the skeleton is 
the union of the solution with a locally optimal solution. In our problem, we again dehne a skeleton 
graph based on a solution and exploit that its treewidth is 0{Vk). This time, the right choice for 
the skeleton seems to be the Voronoi diagram of the k objects forming the solution. Then we exploit 
the fact that this Voronoi diagram has separator cycles of length 0{y/k) to hnd a suitable way of 
separating the instance into subproblems. 

Given a set V of points in the plane, the Voronoi region of a point p G V consists of those points 
of the plane that are closer to p than any other member of V. The boundaries of the Voronoi regions 
are segments that are equidistant to two points from V, forming a diagram that can be considered 
a planar graph. We can define Voronoi regions of a graph and a set D of disjoint connected vertex 
sets in a similar way, by classifying vertices according to the closest object in P. If the graph is 
planar, then we can use the edges on the boundary of the regions (in the dual graph) to construct 
a planar graph that is an analog of the Voronoi diagram. If \T>\ = k, then this diagram has 0{k) 
vertices and, by a well-known property of planar graphs, has treewidth 0{Vk). 

The separator cycle that we need is actually a noose: a closed curve that intersects the graph 
only in its vertices and visits each face at most once. It can be deduced from known results on 
sphere cut decompositions that a planar graph with k vertices has a noose that visits 0(\/k) vertices 
and there are at most |A: faces strictly inside/outside the noose. The noose can be described by a 
cyclic sequence of 0(\/k) vertices and faces. Such a sequence of vertices and faces of the Voronoi 
diagram can be translated into a sequence of shortest paths connecting points from D and the 
branch vertices of the Voronoi diagram, forming a closed cycle in the original graph, separating 
the inside and the outside. These separator cycles are the most important conceptual objects for 
our algorithm. The crucial observation is that, because of the properties of the Voronoi diagram, 
objects of the optimum solution inside the cycle cannot interact with outside world: they cannot 
intersect the cycle and, in covering problems, if a point is sufficiently close to an object inside the 
cycle, then it is sufficiently close also to an object on the boundary. 

Of course, we do not know the Voronoi diagram of the k objects in the solution. But since the 
separator cycles are defined by a selection of 0{y/k) objects/branch vertices, we can enumerate 
^0{Vk) candidates for them. We branch into subproblems indexed by these candidates, and 

in each subproblem we assume that the selected candidate corresponds to a balanced noose in the 
Voronoi diagram of the solution. This assumption has certain consequences and we modify the 
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instance accordingly. For example, we can deduce that certain members of T) cannot be in the 
solution, and hence they can be removed from T). After these modifications, we can observe that the 
problem falls apart into into subproblems: this is because there cannot be any interaction between 
the objects inside and outside the separator. Therefore, we can recursively solve these subproblems, 
where the parameter value is at most |/c. Because of guessing the separator cycle, we eventually 
solve subproblems, each with parameter value at most |fc, which results in the claimed 

^ 0 (\/k) j-m^ning time following from the recursive formula. 

Lower bounds on how k has to appear in the exponent, such as the lower bounds in Theorems|1.9f- 


1.11, can be obtained by parameterized reductions from a W[l]-hard problem. The strength of 


the lower bound depends on how the parameter k changes in the reduction. The reductions based 
on Grid Tiling involve a quadratic blowup and hence are able to rule out o nly algorithms of 
the form f{k)n°^'^^ (assuming ETH). The lower bounds in Theorems 1.9-1.11 are stronger: they 


show that even time algorithms are far from being possible. Therefore, they are very 

different from typical hardness proofs for planar and geometric problems based on Grid Tiling. Of 
particular interest is Theorem [Mil where a tight reduction from Glique seems problematic, as it 
seems difficult to implement the 0{k'^) pairwise interaction of the Glique problem with squares of 
almost the same size. Instead, we rely on a nontrivial hardness result for Partitioned Subgraph 
Isomorphism [3H], which gives a strong lower bound even when the graph H to be found is sparse. 


2 Geometric problems 

Our main algorithmic result is a technique for solving a general facility location problem on planar 
graphs in time With simple reductions, we can use these algorithms to solve 2-dimensional 

geometric problems (see Section 3.2). However, our main algorithmic ideas can be implemented 


also directly in the geometric setting, giving self-contained algorithms for geometric problems. 
These geometric algorithms avoid some of the technical complications that arise in the planar 
graph counterparts, such as the Voronoi diagram having bridges or shortest paths sharing subpaths. 
Unfortunately, a large part of the paper is devoted to the formal handling of these issues. Therefore, 
it could be instructive for the reader to see first a self-contained presentations of some of the 
geometric results. 

Packing unit disks. We start with the Independent Set problem for unit disks: given a 
set V of closed disks of unit radius in the plane, the task is to select a set of k pairwise disjoint 
disks. This problem is known to be solvable in time 13110]. We present another 

algorithm for the problem, demonstrating how we can solve the problem recursively by focusing 
on the Voronoi diagram of a hypothetical solution. This idea appeared recently in the context of 
constructing quasi-polynomial time approximation schemes (QPTAS) for geometric problems (see, 
e.g., I28|), but it has not been used explicitely for exact algorithms. 

While the algorithm we present in this section on its own does not deliver any new result yet, it 
is significantly different from the previous algorithms, which crucially use the notion of “area,” in 
particular, by using the fact that a region of area 0(1) can contain only 0(1) independent unit 
disks. Our algorithm uses only the notion of distance, making it possible to translate it to the 
language of planar graphs, where the notion of area does not make sense. Furthermore, as we 
shall see, generalizations to disks of different radii and to covering problems are relatively easy 
for our algorithm. In what follows, it is oftem more convenient to think in terms of an equivalent 
formulation of the problem where instead of packing disks from V we are packing their centers 
subject to a constraint that every pair of packed centers has to be at distance more than 2 from 
each other. We will switch between these two formulations implicitly. 











(b) Voronoi diagram of 17 points with 3 guard points 
added. Observe that the only infinite regions are 
those of the guard points. 


Figure 1: Voronoi diagrams 


Let "P be a set of points in the plane. The Voronoi region of p G P is the set of those points x in 
the plane that are “closest” to p in the sense that the distance of x and P is exactly the distance of 
X and p (see Figure [^a)). The Voronoi region of p can be obtained as the intersection of half-planes: 
for every p' G P different from p, the Voronoi region is contained in the half-plane of points whose 
distance from p is not greater than the distance from p' . This implies that every Voronoi region is 
convex. 

Even though defining Voronoi diagrams and working with them is much simpler in the plane 
than for their analogs in planar graphs (see Section 4.3), there is a technical difficulty specific to the 
plane. The issue is that the Voronoi region of a point p G P can be infinite and consequently the 
Voronoi diagram consists of finite segments and infinite rays. Therefore, it is not clear what we 
mean by the graph of the Voronoi diagram. While this complication does not give any conceptual 
difficulty in the algorithm, we need to address it formally. 

First, we introduce three new “guard” points into P, at distance more than r from the other 
points and from each other. We introduce these three guards in such a way that every original 
point is inside the triangle formed by them (see Figure [^b)). It is easy to see now that the Voronoi 
region of every original point is finite. Therefore, the only infinite regions are the regions of the 
three guards and it follows that finite segments of the Voronoi diagram form a 2-connected planar 
graph and there are three infinite rays in the infinite face. In the case of the packing problem, if we 
introduce three new disks corresponding to the three guards, then these three disks can be always 
selected into every solution. Thus instead of trying to find k independent disks in the original set, 
we can equivalently try to hnd k + 3 disks in the new set. As increasing /c by 3 does not change the 
asymptotic running time we are aiming for, in the following we assume that the set P of 

center points has this form, that is, contains three guard points. 

If a vertex of the Voronoi diagram has degree more than 3, then this means that there are four 
points appearing on a common circle. In order to simplify the presentation, we may introduce small 
perturbation to the coordinates to ensure that this does not happen for any four points. Moreover, 
we may identify the “endpoints” of the three infinite rays into a new point at infinity. Therefore, in 


9 











the following we assume that the Voronoi diagram is actually a 2-connected 3-regular planar graph. 
Let us emphasize again that these technicalities appear only in the geometric setting and will not 
present a problem when we are proving the main result for planar graphs. 

We are now ready to explain the main combinatorial idea behind the algorithm for finding k 
independent unit disks. Consider now a hypothetical solution consisting of k independent disks and 
let us consider the Voronoi diagram of the centers of these k points (see Figure [^a)). To emphasize 
that we consider the Voronoi diagram of the centers of the k disks in the solution and not the centers 
of the n disks in the input, we call this diagram the solution Voronoi diagram. As we discussed 
above, we may assume that the solution Voronoi diagram is a 2-connected 3-regular planar graph. 
There are various separator theorems in the literature showing that a fc-vertex planar graph has 
balanced separators of size 0{Vk). Certain technicalities appear in these theorems: for example, 
they may require the graph to be triangular or 2-connected, the separator may be a cycle in the 


primal graph or in the dual graph, etc. Therefore, in Section 4.7 we give a short proof showing that 
the known results on sphere cut decompositions imply a separator theorem of the form we need. A 
noose of a plane graph G is a closed curve 6 on the sphere such that 5 alternately travels through 
faces of G and vertices of G and every vertex and face of G is visited at most once. We show that 
every 3-regular planar graph G with k faces has a noose 6 of length 0{y/k) (that is, going through 
0{^/k) faces and vertices) that is face balanced in the sense that there are at most ‘^k faces of G 
strictly inside 5 and at most ifc faces of G strictly outside 5. 

Consider a face-balanced noose 5 of length 0{Vk) as above (see the green curve in Figure ||;b)). 
Noose 5 goes through 0{y/k) faces of the solution Voronoi diagram, which corresponds to a set Q 
of 0{y/k) disks of the solution. The noose can be turned into a polygon F with 0{y/k) vertices the 
following way (see the green polygon in Figure [^c)). Consider a subcurve of 5 that is contained in 
the face corresponding to disk p G Q and its endpoints are vertices x and y of the solution Voronoi 
diagram. Then we can “straighten” this subcurve by replacing it with a straight line segment 
connecting x and the center of p, and a straight line segment connecting the center of p and y. 
Therefore, the vertices of the polygon F are center points of disks in Q and vertices of the solution 
Voronoi diagram. Observe that F intersects the Voronoi regions of the points in Q only. This follows 
from the convexity of the Voronoi regions: the segment between the center oi p G Q and a point on 
the boundary of the region of p is fully contained in the region of p. In particular, this means that F 
does not intersect any disk other than those in Q. 

The main idea is to use this polygon F to separate the problem into two subproblems. Of course, 
we do not know the solution Voronoi diagram and hence we have no way of computing from it the 
balanced noose 5 and the polygon F. However, we can efficiently list candidate polygons. 

By definition, every vertex of the polygon F is either the center of a disk in P or a vertex of the 
solution Voronoi diagram. Every vertex of the solution Voronoi diagram is equidistant from the 
the centers of three disks in T) and for any three such center points (in general position) there is 
a unique point in the plane equidistant from them. Thus every vertex of the polygon F is either 
a center of a disk in T) or can be described by a triple of disks in T). This means that F can be 
described by an 0(\/fc)-tuple of disks from D. That is, by branching into directions, we may 

assume that we have correctly guessed the subset Q of the solution and the polygon F. 

What do we do with the set Q and the polygon F? First, if Q is indeed part of the solution, 
then we may remove these disks from T) and decrease the target number of disks to be found by \Q\. 
Second, we perform the following cleaning steps: 

(1) Remove any disk that intersects a disk in Q. 

(2) Remove any disk that intersects F. 

Indeed, if Q is part of the solution, then no other disk intersecting Q can be part of the solution. 
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(a) The Voronoi diagram of the centers of the disks 
in solution. 



(b) A noose in the Voronoi diagram 




(c) The polygon corresponding to the noose. 


(d) Removing every disk intersecting the polygon 
breaks the problem into two independent parts. 


Figure 2: Using a noose in the Voronoi diagram for divide and conquer. 


Moreover, we have observed above that in the solution the polygon F is contained in the Voronoi 
regions of the points in Q and hence no disk other than the disks in Q intersects F, justifying the 
removal of such disks. We say that these removed disks are banned by ((5,F). 

After these cleaning steps, the instance falls apart into two independent parts: each remaining 
disk is either strictly inside F or strictly outside F (see Figure [^d)). Moreover, recall that the 
noose S was face balanced and hence there are at most |A: faces of the solution Voronoi diagram 
inside/outside 6. This implies that the solution contains at most |A: center points inside/outside F. 
Therefore, in the two recursive calls, we need to look for at most that many independent disks. For 
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k' := 1,..., [^k\, we recursively try to find exactly k' independent disks from the input restricted 
to the inside/outside F, resulting in 2 • recursive callsj^ Taking into account the guesses 

for Q and F, the number of subproblems we need to solve is 2 • |A: • = ^ 0 (\/k) 

otherwise there is no solution) and the parameter value is at most |A: in each subproblem. Therefore, 
if we denote by T(n, k) the time needed to solve the problem with at most n points and parameter 
value at most k, we arrive to the recursion 

r(n, k) = ■ T{n, (2/3)/c). 


Solving the recursion gives 

T(„. k)= ... 

= „o{(i+(|)U(|)2+(§)2+...)v/fc) ^ nO(Vk)^ 

as the coefficient of Vk in the exponent is a constant (being the snm of a geometric series with ratio 
y^2/3). Therefore, the total running time for hnding k independent disks is , This proves the 

first result: packing unit disks in the plane in time Let ns repeat that this result was known 

before [a HQ], but as we shall see, our algorithm based on Voronoi diagrams can be generalized to 
objects of different size, planar graphs, and covering problems. 

Covering points by unit disks. Let ns now consider the following problem: given a set T) of 
unit disks and a set C of client points, we need to select k disks from P that together cover every 
point in C. We show that this problem can be solved in time using an approach based on 

finding separators in the Voronoi diagram. 

Similarly to the way we handled the packing of unit disks, we can consider the Voronoi diagram 
of the center points in the solution. Note, however, that this time the disks in the solution are 
not necessarily disjoint, but this does not change the fact that their center points (which can be 
assumed to be distinct) define a Voronoi diagram. Therefore, it will be convenient to switch to an 
equivalent formulation of the problem described in terms of the centers of the disks: D is a set of 
points and we say that a selected point in T) covers a point in C if their distance is at most 1. 

Similarly to the case of packing, we can try possibilities to guess a set Q C 2? of center 

points and a polygon F that corresponds to a face-balanced noose. The question is how to use F 
to split the problem into two independent snbproblems. The cleaning steps (1) and (2) above for 
the packing problem are no longer applicable: the solution may contain disks intersecting the disks 
with centers in Q and the solution may contain further disks intersecting the polygon F. What we 
do instead is the following. First, if we assnme that Q is part of the solution, then any point in C 
covered by some point in Q can be removed, as it is already covered. Second, we know that in the 
solution Voronoi diagram every point of F belongs to the Voronoi region of some point in Q, hence 
we can remove any point from T) that is inconsistent with this assumption. That is, if there is a 
p € D and n G F such that p is closer to v than to every point in Q, then p can be safely removed 
from D; we say in this case that {Q, F) bans p. For a, p € D and for each segment of F, it is not 
difficult to check if the segment contains such a point v (we omit the details). Thus we have now 
the following two cleaning steps: 

(1) Remove every point from C that is covered by Q. 

^Doing a recursive call for each k' may seem unnecessarily complicated at this point: what we really need is a single 
recursive call returning the maximum number of independent disks, or |fc independent disks, whichever is smaller. 
However, we prefer to present the algorithm in a way similar to how the more general problems will be solved later on. 
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(2) Remove every point from T) that is closer to a point of T than every point in Q. 


Let Din and Dout be the remaining points in V strictly inside/outside L and let Cin and Cout be 
the remaining points in C strictly inside/outside L. We know that the solution contains at most 


|A: center points inside/outside L. Therefore, for k' = 1, 


, [|A:J, we solve two subproblems, with 


point sets (Dm,Cin) and (Dout,Cout)- 

If there is a set of A:in of center points in Din covering Cin and there is a set of /cout center points 
in Dout covering Cout, then, together with Q, they form a solution of \Q\ + fcm + ko^^t center points. 
By solving the defined subproblems optimally, we know the minimum value of kin and fcout required 
to cover Cin and Cout, and hence we can determine the smallest solution that can be put together this 
way. But is it true that we can always put together an optimum solution this way? The problem is 
that, in principle, the solution may contain a center point p G Dout that covers some point q G Cin 
that is not covered by any center point in Din. In this case, in the optimum solution the number of 
center points selected from Din can be strictly less than what is needed to cover Cin and hence the 
way we are putting together a solution cannot result in an optimum solution. 

Fortunately, we can show that this problem never arises, for the following reason. Suppose that 
there is such a p G Dout and q G Cin. As p is outside T and q is inside T, the segment connecting p 
and q has to intersect T at some point u G T, which means dist(p, g) = dist(p, u) + dist(u,g). By 
cleaning step (2), there has to be a p' G Q such that dist(p', v) < dist(p, v), otherwise p would be 
banned and we would have removed it from D. This means that dist(p, q) = dist(p, v) + dist(u, q) > 
dist(p', v) + dist(u, q) > dist(p', q). Therefore, if p covers q, then so does p' G Q. But in this case we 
would have removed q from C in the first cleaning step. Thus we can indeed obtain an optimum 
solution the way we proposed, by solving optimally the defined subproblems. 

As in the case of packing, we have 2 • |A: • subproblems, with parameter value 

at most |A;. Therefore, the same recursion applies to the running time, resulting in an time 

algorithm. 

Packing in planar graphs. How can we translate the geometric ideas explained above to the 
context of planar graphs? Let G be an edge-weighted planar graph and let D be a set of disjoint 
“objects,” where each object is a connected set of vertices in G. Then we can define the analog of 
the Voronoi regions in a straightforward way: for every p G D, let Mp contain every vertex v to 
which p is the closest object in D, that is, dist(u,D) = dist(u,p). By a perturbation of the edge 
weights, we may assume that there are no ties: vertex v cannot be at exactly the same distance 
from two objects pi,P 2 G D. It follows that the sets Mp form a partition of V{G) (here we use 
that the objects in T are disjoint). It is easy to verify that region Mp has the following convexity 
property: if u G Mp and D is a shortest path between v and p, then every vertex of P is in Mp. 

While Voronoi regions are easy to define in graphs, the proper definition of Voronoi diagrams is 
far from obvious and it is also nontrivial how a noose S in the Voronoi diagram defines the analog 
of the polygon T. We leave the discussion of these issues to Section here we only define in an 
abstract way what our goal is and state in Lemma 2.1 below (a simplified version of) the main 


technical tool that is at the core of the algorithm. Note that the statement of Lemma [2.1| involves 
only the notion of Voronoi regions, hence there are no technical issues in interpreting and using 
it. However, in the proof we have to define the analog of the Voronoi diagram for planar graphs 
and address issues such that this diagram is not 2-connected etc. We defer the required technical 
definitions to Section IH 

Let us consider first the packing problem: given an edge-weighted graph G, a set V of d objects 
(connected subsets of vertices), and an integer k, the task is to find a subset D C D of A: pairwise 
disjoint objects. Looking at the algorithm for packing unit disks described above, what we need 
is a suitable guarded separator, which is a pair (Q,r) consisting of a set Q C D of 0{y/k) objects 
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and a subset T C V{G) of vertices. If there is a hypothetical solution C D consisting of k 
disjoint objects, then we would like to have a guarded separator (Q,r) satisfying the following three 
properties: (1) Q is subset of the solution, (2) F is fully contained in the Voronoi regions of the 
objects in Q, and and (3) F separates the objects in in a balanced way. Our main technical result 
is that it is possible to enumerate a set of guarded separators such that for every solution F, 

one of the enumerated guarded separators satisfies these three properties. We state here a simplified 
version that is suitable for packing problems. 

Lemma 2.1. Let G he an n-vertex edge-weighted planar graph, V a set of d connected subsets of 
V{G), and k an integer. We can enumerate (in time polynomial in the size of the output) a set 
Af of d^^'P^ pairs (Q,F) with Q ^V, \Q\ = 0{'/k), F C V{G) such that the following holds. If 
F <LV is a set of k pairwise disjoint objects, then there is a pair {Q,T) G M such that 

1. Q<ZF, 

2. if {Mp)p^jr are the Voronoi regions of F, then F C 

3. for every connected component G of G — T, there are at most objects of F that are fully 
contained in G. 


The proof goes along the same lines as the argument for the geometric setting. After carefully 
dehning the analog of the Voronoi diagram, we can use the planar separator result to obtain a 
noose 5. The same way as this noose was turned into a polygon in the geometric algorithm, we 
“straighten” the noose into a closed walk in the graph connecting using shortest paths 0{y/k) objects 
and 0{^/k) vertices of the Voronoi diagram. The vertices of this walk separate the objects that are 
inside/outside the noose, hence it has the required properties. Thus by trying all sets of 0{y/k) 
objects and 0{pk) vertices of the Voronoi diagram, we can enumerate a suitable set Af. A technical 
difficulty in the proof is that the dehnition of the vertices of the Voronoi diagram is nontrivial. 
Moreover, to achieve the bound instead of we need a nontrivial way of hnding a set 

of dP^^'^ candidate vertices; unlike in the geometric setting, enumerating vertices equidistant from 
three objects is not sufficient. 

Armed with the set N given by Lemma 2.1, the packing problem can be solved in a way analogous 
to how we handled unit disks. We guess a pair (<3,F) G Af that satisfies the three properties of 
Lemma 2.1 From the first property, we know that Q is part of the solution, hence Q can be 


removed from F and the target number k of objects to be found can be decreased by \Q\. From 
the second property, we know that in the solution F the set F is fully contained in the Voronoi 
regions of the objects in Q. Thus we can remove any object from F intersecting F. That is, we 
have to solve the problem on the graph G — T. We can solve the problem independently on the 
components of G — F and the third property of Lemma |2.1| implies that each component contains 
at most |A: objects of the solution. Thus for each component C of G — F containing at least one 
object of F and for k' = 1,..., \_\k\, we recursively solve the problem on G with parameter k', that 
is, we try to find k' disjoint objects in G. Assuming that (Q,F) indeed satisfied the properties 
of Lemma |2.1| for the solution F, the solutions of these subproblems allow us to put together a 
solution for the original problem. As at most d components of G — F can contain objects from F, 
we recursively solve at most d ■ ‘^k subproblems for a given (Q,F). Therefore, the total number of 
subproblems we need to solve is at most d • • |AA| = d ■ ^k ■ (assuming k < d). If 

we denote by T(n, d, k) the running time for |V(G)| < n and \F\ < d, then we arrive to the recursion 
T(n, d, k) = dP'^^'> ■ T{n, d, (2/3)/c) + which gives T{n, k) = d^^P^^ • 

Covering in planar graphs. Let us consider now the following graph-theoretic analog of 
covering points by unit disks: given an edge-weighted planar graph G, two sets of vertices F and C, 
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and integers k and r, the task is to find a set C D of A; vertices that covers every vertex in C. 
Here we say that p £ V covers q £ C if dist(p, q) < r, that is, we can imagine p to represent a ball 
of radius r in the graph with center at p. Note that, unlike in the case of packing, 2? is a set of 
vertices, not a set of connected sets. 

Let be a hypothetical solution. We can construct the set J\f given by Lemma 2^ and guess a 
guarded separator {Q, L) satisfying the three properties. As we assume that Q is part of the solution 
(second property of Lemma 2.1), we decrease the target number k of vertices to select by \Q\ and 
we can remove from C every vertex that is covered by some vertex in Q] let C be the remaining set 
of vertices. Moreover, by the third property of Lemma |2.1t we can assume that in the solution 
the set r is fully contained in the Voronoi regions of the vertices in Q. This means that if there 
is a p £ 'D\Q and u G L such that dist(p, v) < dist(p, Q), then p can be removed from V, since it 
is surely not part of the solution. As in the case of covering with disks, we say that (Q,r) bans 
p. Let T)' be the remaining set of vertices. For every component C of G — S' and A:' = 1,..., [|A:J , 
we recursively solve the problem restricted to C, that is, with the restrictions 22'[C] and C'[C] of 
the vertex sets. It is very important to point out that now (unlike how we presented the packing 
problem above) we do not change the graph G in each call: we use the same graph G with the 
restricted sets 22'[C] and C'[G]. The reason is that restricting to the graph G[G] could potentially 
change the distances between two vertices x,y £ G, as it is possible that the shortest x — y path 
leaves G. 

If kc is the minimum number of vertices in 22'[C] that can cover C'[G], then we know that there 
are \Q\ vertices in 22 that cover every vertex in C. We argue that if there is a solution, we 

can obtain a solution this way. Analogously to the case of covering vertices with disks, we have 
to show that in the solution fF, every vertex in C'\G] is covered by some vertex in 22'[C]: it is not 
possible that there are two distinct components Gi, C 2 of G — T and some q £ C'[Gi] is covered only 
by p G 22'[G 2 ]. Suppose that this is the case, and let P be a shortest p — q path (which has length 
at most r). As p and q are in two different components of G — T, the path P has to intersect T at 
some vertex u G T and we have dist(p, q) = dist(p, v) + dist(u, q). We know that there has to be a 
p' £ Q such that dist(p', v) < dist(p, u), otherwise p would be banned and it would not be in 22'. 
Then the same calculation as in the geometric case shows that p' £ Q also covers q, which means 
that we removed q and it cannot be in C. Thus we can indeed obtain a solution by solving the 
subinstances restricted to the components of G — T. As in the case of packing, we have at most 
d ■ ^k ■ subproblems and the running time follows the same way. 

Covering in planar graphs (maximization version). Let us consider now the variant of 
the previous problem where we want to select k vertices from 22 that cover the maximum number of 
vertices in C. We proceed the same way as before: we guess a guarded separator ((5,r), remove 
from C the set of vertices covered by Q (let m be the number of these vertices), remove from 22 
the set of vertices banned by (QjF), and recursively solve the problem for every component G of 
G — r and every k' = 1,..., L|^J • Having solved these subproblems, we have at our hand the values 
m{G, k'), the maximum number of vertices in C'[G] that can be covered by k' vertices from 22'[G]. 
How can we compute from these values the maximum number of vertices in C that can be covered 
by k vertices from 22? We need to solve a knapsack-type problem: for each component G, we have 
to select a solution containing a certain number 0 < /cc < vertices of 22 such that the sum of the 
/cc’s is at most k — \Q\ and the sum of the values m{G, kc) is maximum possible. Given the values 
m{C, /c'), this maximum can be computed by a standard polynomial-time dynamic programming 
algorithm. Thus instead of just deciding if k vertices from 22 are sufficient to cover every vertex in 
C, we can also find a set of k vertices covering the maximum number of vertices from C. Let us 
point out that in this problem it is really essential that we solve the subinstances for each value of 


15 




k' instead of just finding a minimum/maximum cardinality solution: as the optimum solution may 
not cover all vertices in a component C, computing the minimum number of vertices required to 
cover all vertices in C is clearly not sufficient. 

Covering in planar graphs (nonuniform radius). A natural generalization of the covering 
problem is when every vertex p G D is given a radius r(p) > 0 and now we say that p covers a vertex 
g G C if dist(p, q) < r{p). That is, now the vertices in D represent balls with possibly different radii. 

There are two ways in which we can handle this more general problem. The first is a simple 
graph-theoretic trick. Let R be the maximum of all the r(p)’s. For every p € R, let us attach a 
path of length R — r{p) to p, let us replace p mV with the other end p' of this path. Observe that 
now a vertex g G C is at distance at most r(p) from p if and only if it is at distance at most R 
from p'. Therefore, by solving the resulting covering problem with uniform radius R, we can solve 
the original problem as well. Note that, interestingly, this trick uses the flexibility of the problem 
being stated in terms of graphs and there is no analogous geometric trick to solve the problem 
of covering points with disks of nonuniform radius in the plane: by attaching the paths, we are 
distorting the distance metric in a specific way, whose geometric meaning is difficult to interpret in 
the 2-dimensional plane. 

The second way is somewhat more complicated, but it seems to be the robust mathematical 
solution of the issue and it has also a geometric interpretation. The issue of nonuniform radius can 
be handled by working with the additively weighted version of the Voronoi diagram, that is, instead 
of defining the Voronoi regions of V by comparing the distances dist(p, u) for p €V, we compare 
the weighted distances dist(p, v) — r{p). It can be verified that the main arguments of the algorithm 
described above go through. For example, it remains true that Voronoi regions have the convexity 
property that if v is in the region of p, then the shortest path from u to p is fully contained in the 
Voronoi region of p. Given a guarded separator (QjF), the crucial property allowing us to separate 
the problem to the components of G — F was that if p G P and q £ C are in two different components 
and p covers q, then some p' £ Q also covers q. This property also remains true: we can redo the 
same calculation to compare dist(p, q) — r{p) and dist(p', q) — r{p'). 

Combining covering and packing. Comparing the algorithms for packing objects in planar 
graphs and for covering vertices by balls of radius r in planar graphs, we can observe that the set V 
played very different roles in the two algorithms: in the packing algorithm V contained the actual 
objects, while in the covering algorithm V contained only the centers of the balls. The reason for 
this is fundamental: in order to define the Voronoi regions of the solution T" C T>, we need F to be 
a set of disjoint objects (which is not true for the balls in the solution of the covering problem, but 
of course true for the centers of these balls). This means that we cannot generalize the covering 
algorithm to objects different from metric balls simply by putting objects into V more general than 
single vertices. However, we have no trouble obtaining a common generalization if we require that 
the objects selected from V are disjoint. That is, in the independent covering problem, we are given 
a set V of connected objects in a planar graph G, a set C of vertices, and integers k and r, the task 
is to select a pairwise disjoint subset T" C T> of size exactly k that covers the maximum number 
of vertices in C (where, as usual, an object p G V covers a vertex g G C if dist(p, q) < r). Exactly 
the same algorithm as above goes through: as the solution T is disjoint, it is possible to define the 
Voronoi regions. 

In the next section, we define the Disjoint Network Coverage problem that generalizes 
all our applications. The main algorithmic result is expressed as an algorithm for this problem 
(Theorem |3.1[ ). To handle covering with nonuniform radius, the problem formulation allows different 
radius r(p) for each p G V. However, if we allow nonuniform radius, then the requirement stating 
that we have to select disjoint objects should be replaced by a technical condition that we have 
to select a normal family of objects (see Section 3.1). Essentially, this condidition states that in 
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the solution T, every vertex of every object p £ J- should be in the weighted Voronoi region of 
p. In other words, if we allow arbitrary radii, it could in principle happen that one object p has 
so large radius compared to some other object p' that the weighted Voronoi region of p actually 
contains some vertices of p', thus “eating” the object p' itself and making it not contained in its 
own Voronoi region. The normality condition states that, in addition to being disjoint, the solution 
F does not contain snch pathological situations. Observe that when the objects are single vertices 
only, then this problem is not very dangerous: if p “dominates” p' in the way described above, 
then every client covered by p' is also covered by p, and we can discard p' from T> because taking p 
instead is always more prohtable. Unfortunately, if the objects are not just single vertices, or they 
are equipped with nonuniform costs, then this technicality makes the arising independent covering 
problems with nonuniform radii somewhat unnatural. 


3 The general problem 


In this section we introduce the generic problem Disjoint Network Coverage. The main 
result of this paper is an algorithm for this problem, expressed in Theorem 3.1 Before we give the 


algorithm for DISJOINT Network Coverage, in Section [3^ we shall see how the concrete results 
mentioned in the introduction (i.e., Theorem 1.1 1.7) follow from Theorem 3.1 by simple reductions. 


3.1 Problem definition 

Suppose we are given an undirected graph G embedded on a sphere S, together with a positive 
edge weight function w: E{G) —)• M"*". We are given a family of objects T> and a family of clients C. 

Every object p € T> has three attributes. It has its location loc(p), which is a nonempty subset 
of vertices of G such that G[loc(p)] is a connected graph. Moreover, it has its cost \{p), which is 
a real number. Note that costs may be negative. Finally, it has its radius r{p), which is just a 
nonnegative real value denoting the strength of domination imposed by p. Note that locations of 
objects can intersect, and even there can be multiple objects with exactly the same location. 

Every client q (z C has three attributes. It has its placement pla(q), which is just a vertex of 
G where the client resides. It has also its sensitivity s{q), which is a real value that denotes how 
sensitive the client is to domination from objects. Finally, it has also prize 7r(g), which is a real 
value denoting the prize gained by dominating the client. Note that there can be multiple clients 
placed in the same vertex, and the prizes may be negative. 

We say that a subfamily C D is normal if locations of objects from F are disjoint, and moreover 
dist(loc(pi), loc(p 2 )) > r{pi) — r{p 2 ) for all pairs {pi,P 2 ) of different objects in F] here dist(X, U) 
denotes the distance between two vertex sets in G w.r.t. edge weights w. As we require the same 
inequality for the pair {p 2 ,Pi) as well, it follows that actually dist(loc(pi), loc(p 2 )) > |j’(pi) — i"iP 2 )\ 
is true. In particular, this implies disjointness of locations of objects from F, but if all radii are 
equal, then normality boils down to just disjointness of locations. Note also that a subfamily of a 
normal family is also normal. 

We say that a client q is covered by an object p if the following holds: dist(pla(g), loc(p)) < 
s{q) + r{p). In other words, a client gets covered by an object if the distance from the object does 
not exceed the sum of its sensitivity and the object’s domination radius. 

We are finally ready to define Disjoint Network Coverage. As input we get an edge-weighted 
graph G embedded on a sphere, families of objects V and clients C (described using locations, costs, 
radii, placements, sensitivities, and prizes), and a nonnegative integer k. For a subfamily Z <^T>, we 
define its revenue, denoted n(Z), as the total sum of prizes of clients covered by at least one object 
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from Z minus the total sum of costs of objects from Z. In the Disjoint Network Coverage 
problem, the task is to find a subfamily Z ^T) such that the following holds: 

(i) Family Z is normal and has cardinality exactly k. 

(ii) Subject to the previous constraint, family Z maximizes the revenue Il{Z). 

It can happen that there is no subfamily Z satisfying property |(i)[ In this case, value —oo should 
be reported by the algorithm. For an instance {G,'D,C, k) of Disjoint Network Coverage, we 
denote d = \V\, c = \C\, and n = |D(G)|. Observe that we can assume that k < d and that the 
input graph G is simple: loops can be safely removed, and it is safe to only keep the edge with the 
smallest weight from any pack of parallel edges. 

Note that if we do not provide any clients in the problem and we set all the radii to be equal 
to 0, then we arrive exactly at the problem of packing disjoint objects in the graph. However, by 
introducing also clients we can ask for (partial) domination-type constraints in the graph, as well as 
define prize-collecting objectives. In this manner. Disjoint Network Coverage generalizes both 
packing problems as well as covering problems. The caveat is, however, that we need to require that 
the objects that cover clients are pairwise disjoint, and moreover they have to form a normal family. 

The main result of this paper is the following theorem. 

Theorem 3.1 (Main result). Disjoint Network Coverage can be solved in time dP^'^^ • 


3.2 Applications of Theorem 3.1 


In this section we provide formal argumentation of how the concrete results mentioned in the 
introduction follow from Theorem 3.1 While for problems on planar graphs the reductions are 


straightforward, for geometric problems some non-trivial technicalities arise. We remark that we 
are proving the exact statements given in the introduction, but the generality of the Disjoint 
Network Coverage problem would equally easily allow us to solve more general variants with 
different costs, prizes, sensitivities of clients, etc. 


Applications to planar graphs: Theorems OB 

Restated Theorem 1.1 (packing connected sets). Let G be a planar graph, T) be a family of 
connected vertex sets of G, and k be an integer. In time • n^^^l, we can find a set of k 

pairwise disjoint objects in V, if such a set exists. 

Theorem o follows by taking V to be the set of objects (identified with their locations), 
assigning each of them radius and cost equal to 0, and considering no clients (putting C = 0). As 
distances play no role in the problem, we can simply put unit weights. 


Restated Theorem 1.2 (covering vertices with centers of different radii). Let G be a 

planar graph, let D,G Q V{G) be two subset of vertices, let r: D —)• Z"*“ be a function, and k be an 
integer. In time ■ n^^^\ we can find a set S D of k vertices that maximizes the number 

of vertices covered in C, where a vertex u £ G is covered by v G S if the distance of u and v is at 
most r(v). 


For Theorem 1.2, for every v £ D we construct an object py with loc(p^) = {u}, A(u) = 0 


and r{py) = r{v). Then, for every u £ G we construct a client with pla(g„) = u, s{qu) = 0 
and '7T{qu) = 1- Let V,C be the sets of constructed objects and clients, respectively. We claim 


18 







that the optimum value for the input instance is equal to the maximum among optimum revenues 
for instances {G,T>,C,k') of Disjoint Network Coverage, for 0 < k' < k. On one hand, any 
solution to any such instance {G,'D,C,k') trivially yields a solution to the input instance with the 
same revenue. On the other hand, if 5 C D is a solution to the input instance, then observe that 
without loss of generality we may assume that S does not contain any two distinct vertices vi,V 2 
such that dist(ui,i; 2 ) < r{vi) — r{v 2 )- This is because when dist(ui,U 2 ) < t(ui) — r{v 2 ), then any 
client covered by vi is also covered by V 2 , and vi can be safely removed from S. Then a solution 
S with this property corresponds to a subfamily Z ^T) that is normal and has size \S\ < k, i.e., 
it is a solution for {G,'D,C, |5|) with revenue equal to the number of vertices of C covered by S. 
We conclude that to solve the input instance it suffices to solve all the instances {G,V,C, k') of 
Disjoint Network Coverage for 0 < k' < k, using the algorithm of Theorem |3.l[ 


Restated Theorem 1.3 (covering vertices with independent objects). Let G be a planar 
graph, let V be a set of connected vertex sets in G, let G (G) be a set of vertices, and let k be 
an integer. In time • n^^^\ we can find a set S of at most k pairwise disjoint objects in D 

that maximizes the number of vertices of G in the union of the vertex sets in S. 


Theorem 1.3 follows by taking D to be the set of objects (identified with their locations), 
assigning each of them radius and cost equal to 0, and for every u G G introducing a client Qu with 
pla((jrM) = u, s{qu) = 0 and 7r{qu) = 1. As distances play no role in the problem, we can simply put 
unit weights. 


Applications to geometric problems: Theorems |1.4t]1.7| . In all the geometric reductions 
that follow, we perform arithmetic operations only of the following types: 

• Computing new points on the plane with coordinates given as constant-size rational expressions 
(i.e., yielded by the four basic arithmetic operations) over the input coordinates and radii. 

• Computing Euclidean distances between the obtained points. 

Note that finding intersection of two segments between pairs of points given on the input gives a 
point that is compliant with this requirement. 

Thus, if the input coordinates and radii are s-bit integers, then the coordinates of the points 
obtained in the reductions can be stored as rationals represented using poly(s) bits. Then the 
distances can be safely stored using poly(s, n)-bit precision in order to avoid rounding errors when 
adding at most n of them. Similarly, if the input coordinates and radii are given as floating point 
numbers with s bits before the point and s bits after the point, then we can scale them up to get 
2s-bit integers, and perform the same analysis. 

Restated Theorem 1.4 (packing disks). Given a setT> of disks (of possibly different radii) in 
the plane, in time we can find a set of k pairwise disjoint disks, if such a set exists. 

Proof. We provide a reduction to the problem considered in Theorem o Let I = {V, k) be the 
input instance. Based on the disk set V, we define a planar graph G as follows. First, construct a 
set of points X (see Figure]^ by taking (a) all the centers of disks from V, and (b) for any two 
disks Ai,A 2 that intersect, an arbitrarily chosen point in their intersection (for example, a point on 
the segment between the centers for which distances to the centers are in the same ratio as radii 
of the disks). Start the construction of G by putting V{G) = X. Then, for each pair of points 
xi,X 2 G X, draw a segment between xi and X 2 on the plane, provided that no other vertex of X 
lies on this segment. These segments form so far the edge set of G, but they may cross. To make G 
planar, for every crossing of two segments introduce a new vertex at the intersection point, and 
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Figure 3: Proof of Theorem 1.4 The black dots are the centers of the disks, the red dots 


the points chosen in the intersections of pairs of disks, and the blue dots were introduced at 
intersections of the segments to make the graph planar. 


are 

the 


connect all the old and new points according to the segments. That is, every former edge between 
two points xi^X 2 ^ X becomes a path from xi to X 2 traversing consecutive crossing points, and the 
embedding of this path is exactly the segment between xi and X 2 - To make G edge-weighted, to 
every edge uv G E{G) assign weight equal to the Euclidean length of the segment between u and v. 
Observe that G is connected and |P(G)| < 0{\V\‘^). 

Now, for every disk A € D with radius construct a vertex set X{A) as follows: X{A) comprises 
all the vertices of G that are at distance at most from the center of A (which belongs to 14(0)) 
in the graph G. Obviously G[X{A)\ is connected, because it is defined as a ball in the graph G. 
Moreover, by the triangle inequality it follows that every vertex of X{A) is actually embedded 
into the disk A. Note, however, that X(A) may not contain some intersection points that are 
actually embedded into the disk A, but in G they are at distance larger than from the center of 
A. Nonetheless, X[A) contains all the elements of X that are embedded into the disk A. This is 
because by the construction of G, for any x G X n ^ we have that x is connected to the center of A 
via a path consisting of parallel segments, which in particular has length equal to the Euclidean 
distance from x to this center. 

Let Vq = {X{A) I A G V} be the family of constructed vertex subsets, and let I' = {G, Vq, k) 
be the constructed instance of the problem considered in Theorem 0 We now claim that the input 
instance I has a solution if and only if the constructed instance I' has a solution. 

On one hand, if 5 C T> is a subfamily of k disjoint disks, then Sq = {X{A) \ ^ G 5} is a 
subfamily of k vertex sets that are pairwise disjoint, and hence I' has some solution. On the other 
hand, suppose there exists some solution Sq to and let 5 C P be the corresponding set of k 
disks. We claim that these disks are pairwise disjoint. Suppose that, on the contrary, there are 
some distinct disks Ai,A 2 G 5 that intersect. Recall that for these disks we have added some point 
X G Ai n A 2 to X. As we have argued, by the construction of G this means that x belongs both 
to X{Ai) and to X(A 2 ), which is a contradiction with the fact the vertex sets in Sq are pairwise 
disjoint. □ 
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Figure 4: Proof of Theorem 1.5 The black dots are the vertices of the polygons and the red dots 
are the points introduced at the intersection of segments. The red segments were introduced to 
make the graph connected. 


Restated Theorem 1.5 (packing simple polygons). Given a sefD of simple polygons in the 
plane, in time we can find a set of k polygons in D with pairwise disjoint closed 

interiors, if such a set exists. Here n is the total number of vertices of the polygons in V. 

Proof. We give a reduction to the problem considered in Theorem o Let I = {V, k) be the input 
instance. Based on the set of polygons T), define a planar graph G as follows. Start with taking 
V{G) to be the set of all the vertices of all the polygons of T), and then for every polygon A € P 
draw all the sides of A as edges of G, embedded as respective segments (see Figure Q. Then perform 
a similar construction as in the proof of Theorem 1.4 in order to make G planar, for every crossing 
of two segments introduce a vertex on their intersection. The edge set of G is defined as the set 
of all the segments xy between two vertices x,y € V{G), being either the original vertices of the 
polygons or the introduced crossing points, such that xy is contained in some side of an original 
polygon and the interior of xy does not contain any other vertex of V{G). Thus, this definition 
yields a planar embedding of G. Note that the definition is also correct when we have two sides xy 
and x'y' of two different polygons that share a common subinterval (not being a single point); then 
vertices x, y count as crossing points subdividing x'y' and vice versa. 

Observe that G defined in this way is planar, but may not be connected. Hence, we make 
G connected by iteratively taking two connected components Ci, C 2 of G that are incident to a 
common face, and adding an edge within this face between two arbitrarily chosen vertices from Ci 
and C 2 . Observe that |F(G)| = 0{n^). 

Now, for every polygon A £P construct a vertex set X{A) as follows: X(A) comprises all the 
vertices of G that are embedded into the polygon A. We regard polygons as closed, so X{A) in 
particular contains all the vertices of A. Note that the perimeter of A naturally induces a simple 
cycle in G, and X{A) is exactly the set of points enclosed or on this cycle. Since G is connected, it 
also follows that G[X{A)] is connected. Let Pq = {X{A) \ A G Pj be the set of constructed vertex 
sets, and let I' = {G,Po, k) be the constructed instance of the problem considered in Theorem 0 
We now claim that the input instance I has a solution if and only if the constructed instance I' has 
a solution. 

On one hand, if 5 C is a subfamily of k disjoint polygons, then Sq = {X{A) | H G 5} is 
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a subfamily of k vertex sets that are pairwise disjoint, and hence I' has some solution. On the 
other hand, suppose there exists some solution So to and let 5 C P be the corresponding set of 
polygons. We claim that these polygons are pairwise disjoint. Suppose that, on the contrary, there 
are some distinct polygons Ai,A 2 G S that intersect. If the perimeters of Ai and A 2 intersect, then 
each their intersection point x belongs to V{G) and is contained in X{Ai) n X(A 2 ), contradicting 
the fact that X{Ai) and ^(^ 2 ) are disjoint. Otherwise, either Ai is entirely contained in A 2 , or A 2 
is entirely contained in Ai. In the former case every vertex of Ai belongs to V{G) and is contained 
in X{Ai) n X(A 2 ), and in the latter case every vertex of A 2 has this property. In both cases we 
obtain a contradiction with the fact that X{Ai) and X{A 2 ) are disjoint. □ 

Restated Theorem 1.6 (covering points with disks). Given a set C of points and a set D of 
disks (of possibly different radii) in the plane, in time we can find a set of k disks 

in V maximizing the total number of points they cover in C. 

Restated Theorem 1.7 (covering points with sqnares). Given a set C of points and a set D 

of axis-parallel squares (of possibly different size) in the plane, in time we can find 

a set of k squares in T> maximizing the total number of points they cover in C. 


Proof. We prove both theorems at the same time by showing that the algorithm can be constructed 
whenever sets from V are defined as balls in some norm 11 • 11 on . More precisely, there exists a 


norm 


on 


% such that every A G V is defined as = {x G | \\x — c(A)|| < r{A)}, where 
c{A) is the center of ball A and r[A) is its radius. Then Theorem 1.6 follows by taking || • || to be 
the ^ 2 -iiorm and Theorem 1.7 follows by taking || • || to be the ^oo-norm. 

We give a reduction to the problem considered in Theorem |1.2| and to this end we perform a 


similar construction as in the proof of Theorem 1.4 Let I = (V,C,k) be the input instance. Let X 
be the set of all the centers of balls from P and all the points from C. Start the construction of G by 
putting V{G) = X. Then, for each pair of points xi,X 2 G X such that xi is a center of a ball from 
P and X 2 is a point from C, draw a segment between xi and X 2 on the plane, provided that no other 
vertex of X lies on this segment. These segments form so far the edge set of G, but they may cross. 
To make G planar, for every crossing of two segments introduce a new vertex at the intersection 
point, and connect all the old and new points naturally according to the segments. To make G 
edge-weighted, observe that every edge e of G is embedded as a segment on the plane, so put w(e) 
to be the || • ||-length of this segment. Observe that G is connected and |1^(G)| = 0{\P\‘^\C\‘^). 

Now, let us set D = {c{A) \ A G P} and G = C. Construct an instance I' of the problem 
considered in Theorem |1.2| by taking graph G with subsets of vertices D and C, assigning r{c{A)) = 
r{A) for each c{A) G D, and putting budget k. We now claim that a vertex x G C is covered by some 
c{A) G D in the instance I' if and only if point x G C belongs to ^ G 2? in the input instance I. To 
prove this it suffices to show that distG(®, c{A)) = ||x — c(^)||, since x is covered by c{A) in I' if and 
only if distG(a;, c(A)) < r{A), whereas x belongs to ^ in / if and only if ||x — c(74)|| < r{A). First, 
by the triangle inequality it follows that distances in G are lower bounded by distances w.r.t. norm 
II • II on the plane, so distG(x,c(^)) > ||x —c(A)||. On the other hand, since x,c{A) G X, then by the 
construction of G there exists a path in G that connects x and c{A) and consists of parallel segments. 
In particular, this path has total length ||x — c(^)||, which means that distG(3:,c(^)) < ||x — c(74)||. 
We conclude that indeed distG(a;, c(^)) = ||x — c(^)||. 

Solutions to the input instance I correspond one-to-one to solutions of the constructed instance 
I', and the argumentation of the previous paragraph shows that this correspondence preserves the 

on 

□ 


number of points covered by the solution. Hence, it suffices to run the algorithm of Theorem |1.2 
the instance I'. 


22 







Figure 5: Proof of Theorem 1.6 The black dots are the centers of the disks in T), the red 


in the points in C, and the blue dots were introduced at the intersections of the segments 
the graph planar. 


dots are 
to make 


4 The main algorithm: proof of Theorem 3.1 


4.1 Notation and general definitions 


We first establish notation and recall some general definitions that will be used throughout the 
proof. 

For a graph G, by V{G) we denote its vertex set and by E{G) we denote its edge set. Graph G 
is simple if every edge connects two different vertices, and every pair of vertices is connected by at 
most one edge. On the other hand, G is a multigraph if we allow (a) parallel edges, that is, multiple 
edges connecting the same pair of vertices, and (b) loops, that is, edges that connect some vertex u 
with itself. A graph induced by a vertex subset X C V{G), denoted G[A], has X as the vertex set, 
and its edge set comprises all the edges of G whose both endpoints belong to X. A graph spanned 
by an edge subset F C E{G), denoted G[F], has F an the edge set, and its vertex set comprises all 
the vertices of G incident with at least one edge of F. 

In this work we will be mostly working with edge-weighted graphs. That is, we assume that the 
input graph G is given together with a positive edge weight function w: E{G) —)> M"*". The distance 
between two vertices u,v £ V{G), denoted dist(M, u), is defined as the minimum total weight of a 
path connecting u and v. We extend this notation to subsets in the obvious manner: dist(M, T) is 
the minimum weight of a path from u to any vertex of Y, and dist(X, T) is the minimum weight 
of a path from any vertex of X to any vertex of Y. Note that the dist(-, •) function satishes the 
triangle inequality: dist(u,u) -|-dist(u, rc) > dist(u,u;), and similar inequalities hold for sets as well. 

Let S be the sphere = {(x, y, z) £ \ + z'^ = 1}. A curve on S is a homeomorphic 

image of the interval [0,1], and a closed curve on S is a homeomorphic image of the circle 
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= {{x,y) G I = 1 }. Note that, in particular, a curve and a closed curve do not have 

self-crossings, i.e., no point of S is visited twice. The Jordan curve theorem states that for any 
closed curve 7 on S, E \ 7 consists of two connected sets homeomorphic to open disks. A sphere 
embedding of a graph G is a mapping that maps vertices of G to distinct points of E, and edges 
of G to curves connecting respective endpoints that do not intersect apart from their endpoints. 
In case e G E{G) is a loop, the corresponding curve connects a point with itself, i.e., it is a closed 
curve. Note that this definition is also valid for multigraphs. It is well known that a graph is planar, 
i.e., it can be embedded on the 2 -dimensional plane in the same manner, if and only if it can be 
embedded on a sphere. 

If G is a sphere-embedded (multi)graph, then a face of G is an inclusion-wise maximal connected 
set that does not contain any point of the embedding of G. Note that in particular every face is 
open. It is well known that if G is connected, then every face is homeomorphic to an open disk. 
Note, however, that if G has bridges or cutvertices, then there may be faces whose closures are not 
homeomorphic to closed disks. If the closure of each face is homeomorphic to a closed disk, then we 
say that the embedding is a 2-cell embedding. The set of faces of a (sphere-embedded) graph G 
is denoted by F{G). The well-known Euler formula states that for every connected multigraph G 
embedded on a sphere, it holds that |E(G)| — |E(G)| -|- |T(G)| = 2 . 

Let / G F{G) be a face of a connected, sphere-embedded multigraph G. We say that a vertex v 
or an edge e is incident with / if it is contained in the closure of /. By the boundary of /, denoted 
df, we mean a walk in G that visits consecutive vertices and edges incident with / in the order of 
their appearance around /. Note that if the closure of / is homeomorphic to a closed disk, then 
df is a simple cycle in G. However, otherwise a vertex can be visited more than once on df (this 
can happen if it is a cutvertex in G), and an edge can be traversed more than once on df (this 
can happen if it is a bridge in G). In case u is a cutvertex of G, then there is a face / of G that 
appears multiple times around v in the planar embedding. Whenever some curve is leaving v to f 
(or entering v from /), we will distinguish these different appearances and call then directions. Note 
that every such direction corresponds to an appearance of v on df. 

A sphere-embedded graph is triangulated, if for every face /, walk df is a triangle in G; in 
particular, the graph is simple and the embedding is a 2-cell embedding. It is well-known that for 
any simple graph G embedded on a sphere, one can add edges to the embedding so that the new 
supergraph is triangulated. 

4.2 Simplifying assumptions 

Throughout the proof we assume that {G,V,C,k) is the input instance of Disjoint Network 
Coverage. Before we proceed to the proof, let us make a few simplifying assumptions about the 
input instance {G,V,C,k). These assumptions can be made without loss of generality, and they 
will streamline further reasonings. 

Firstly, we will assume that the graph G is connected. Indeed, for disconnected graph we 
can apply the algorithm for every connected component of G separately, for all the parameters 
i = 0 , 1 , 2 ,...,/?, and then merge the results using a simple dynamic programming algorithm. 

Secondly, we will assume that all the values dist(u, v) and dist(u, v) — r{p) for u,v £ V (G) and 
p £ T> are pairwise different, and moreover that the shortest paths between pairs of vertices in G are 
unique. This can be easily obtained using the standard technique of breaking ties by adding very 
small, distinct values to the weights of edges as well as to the radii of objects. To do this, we need 
to use {d -|- more bits of precision in the representation of floating point numbers. Since we 

never estimate precisely the exact running time of polynomial-time subroutines run on G, this does 
not influence the claimed asymptotic running time of the algorithm. For brevity, we will denote this 
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assumption by (ijk). 

Thirdly, we fix some sphere embedding of G, and without loss of generality we assume that G 
is triangulated. This can be achieved by triangulating every face of G using edges of weight +oo 
(equivalently, a value much larger than all the radii and sensitivities). After this operation, every 
face of G is a triangle whose closure is homeomorphic to a closed disc (i.e., it is a 2-cell embedding). 
Note that thus G is 2-connected, so in particular it does not have any cut-vertices, bridges, or 
vertices of degree 1. We also assume that G has more than 3 vertices, since we can always add a 
new vertex within a triangular face connected with its vertices via edges of weight -|-oo. 

Finally, for every object p £ D we arbitrarily distinguish one vertex cen(p) G loc(p), called the 
center oi object p. We also pick an arbitrary spanning tree T{p) of G[loc(p)], rooted at cen(p). For 
concreteness, the reader may think of T{p) as the tree of shortest paths from cen(p) in G[loc(p)], 
but we will never use this property. 


4.3 Voronoi partitions 

We want to dehne analogs of Voronoi regions and related notions in the context of planar graphs; 
given a set J- of objects, we would like to partition the vertices of the graph according to which 
object in T is closest to a vertex. For this to make sense, we clearly need the assumption that the 
objects in T are pairwise disjoint, otherwise it is not clear how to classify vertices in the intersection 
of the objects. Moreover, to express the fact that different objects have different power of domination 
(r(p) can be different for different objects), we need additively weighted Voronoi regions. That 
is, instead of comparing the distances to locations of objects from we compare the values of 
distances to loc(p) minus r{p), for p £ T. Here we need the technical assumption that the family T 
is normal, which rules out degenerate situations by ensuring that every point of the object is in the 
Voronoi region of the object itself, and hence the Voronoi regions are nonempty and connected. 

Formally, suppose C P is a normal subfamily of objects. Then we define the Voronoi partition 
of G with respect to T, denoted Mj-, as a partition of V(G) into parts {Mpjpgj- as follows: a vertex 
V £ V{G) belongs to part Mp^ if and only if dist(u, loc(po)) — fivo) is the smallest value among 
{dist(u, loc(p)) — r(p)}pgj-. Note that assumption (ijk) implies there can be no ties between different 
objects. Set Mp shall be called the Voronoi region of p. For some Voronoi partition Mj- and p £ iF, 
the Voronoi region of p will be also denoted by Mj-(p). Whenever the family we are working with 
will be clear from the context, we shall drop the index F. This remark holds for all the objects 
defined in the sequel that depend on the currently considered family F. 

The following simple lemma shows that the normality condition is a sanity check that the 
Voronoi regions are well-defined. 

Lemma 4.1. Suppose F <£T) is a normal subfamily of objects and let M = Mj- be the Voronoi 
partition w.r.t. F. Then for every p £ F, it holds that loc(p) C M.{p). 

Proof. For the sake of contradiction, suppose there are two distinct objects pi and p 2 such that 
some vertex v £ loc(pi) belongs to M(p 2 )- Then in particular we have that 

dist(u, loc(p 2 )) -r{p 2 ) < dist(u, loc(pi)) -r{pi) = -r{pi), 


and hence 

dist(loc(pi), loc(p 2 )) < dist(u, loc(p 2 )) < r{p 2 ) — r{pi). 

This is a contradiction with the definition of normality. □ 

The following simple claim shows that Voronoi regions behave as expected. 


25 



Lemma 4.2. Suppose T ^T) is a normal subfamily of objects and let M = Mj- he the Voronoi 
partition w.r.t. T. Then for every p € J- and every v G the unique shortest path from v to 

loc(p) is entirely contained in M(j?). In particular, G[M(p)] is connected. 

Proof. Let P be the shortest path connecting v with some w G loc(p). Note that the uniqueness 
of P follows from (Jjk). For the sake of contradiction, suppose there exists v' G V{P) such that 
dist(?;', loc(p')) — r{p') < dist(n', loc(p)) — r{p), for some p' £ if, p' ^ p. Since P is the shortest 
path, we have: 


dist(t', loc(p)) 


r{p) = dist(r;, w) — r{p) 

= dist(?;, ?;') + dist(n', rc) — r(p) 

= dist(r;, v) + dist(n', loc(p)) — r{p) 

> dist(r;, r;') + dist(n', loc(p')) — r(p') > dist(t>, loc(p')) — r(p'). 


This is a contradiction with v G M(p). The second part of the claim follows from the first part. 
Lemma 4.1 and the fact that G[loc(p)] is connected. □ 


Lemma |4.2| gives rise to the following definition of a spanning tree of a Voronoi region. For a 
normal subfamily of objects jF TV and some p £ J-, we shall distinguish a spanning tree Tjr[p) of 
G[M(p)] constructed as follows: 


• start with Tjr{p) = T{p)-, 

• for every v £ M(p), add to Tjrfp) the unique shortest path from v to loc(p). 


By Lemma 4.2 and properties of shortest paths, Tjrfp) dehned in this manner is a spanning tree of 
G[M(p)]. Note that we may view Tjrfp) as rooted in cen(p), the center of p. However, Tjr[p) is not 
necessarily the tree of shortest paths from cen(p) in G[M(p)]; for example, it is possible that the 
shortest path between two vertices of loc(p) is not contained fully in loc(p). 

The following statement is straightforward, but it is instructive to state it explicitly. 


Lemma 4.3. Let V' T IF TV he two normal families of objects. Then for every p £ IF' it holds 
that Mj-/(p) D Mj-(p) and Tjri(p) is a supergraph ofTjr{p). 


In other words, when we restrict the family of objects, then the Voronoi regions and their 
spanning trees can only grow. 


4.4 Voronoi (pre-)diagrams 

Let us fix some normal subfamily of objects IF, let I = IT"!, and let M = Mj- be the corresponding 
Voronoi partition. Intuitively, M partitions G into regions, such that if we put an edge between 
adjacent regions, then the resulting graph is planar. We now formalize this concept by introducing 
Voronoi prediagrams and Voronoi diagrams. More precisely, Voronoi (pre-)diagrams correspond to 
the dual of the adjacency graph between the regions. 

A Voronoi prediagram for T" is a subgraph FLjr of G*, the dual of G, defined as follows: 

1. Start with the graph G*. 

2. For every p £ IF, remove all the dual edges corresponding to edges of T{p). 

3. Iteratively remove vertices of degree 1 in the obtained graph, up to the point when the 
minimum degree of the graph is 2. 

In this definition we drop the weights of edges, i.e., H is an unweighted graph. See Figure for an 
example. The following lemma formalizes the properties of the Voronoi prediagram. 
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Figure 6: Construction of the prediagram in an area between three objects: red, yellow, and green. 
The Voronoi regions are depicted in light red, light yellow, and light green. Gray edges belong to 
the primal graph G, and thick gray edges belong to spanning trees Tjr[p). The white vertices and 
thick blue edges belong to the prediagram T-Lj^. 


Lemma 4.4. The prediagram TL is a simple, connected graph embedded on the sphere, where every 
vertex is of degree 2 or 3. Moreover, TL has exactly I faces, and faces ofTL correspond to objects of 
T in the following sense: object p is associated with a face f* G FiTl) such that f* contains all the 
vertices ofWl(p), and no other vertex. 

Proof. Observe that TL is a, subgraph of G*, so it inherits a sphere embedding from G*. The fact 
that TL is simple follows from the fact that G* is simple, since G is triangulated, simple, and has 
more than 3 vertices. Also, the fact that G is triangulated implies that G* is a 3-regular graph. 
Hence TL, as a subgraph of G*, has maximum degree at most 3. The fact that G* has minimum 
degree at least 2 follows from the construction. 

We now show that TL is connected. Obviously G* is connected, as a dual of a sphere-embedded 
graph. We now argue that during the step when the dual edges to the edges of some T{p) are 
removed, we also cannot spoil the connectivity of TL. Let VF be a walk in G* that traverses the set 
of faces incident to vertices and edges of T{p), but does not use any dual edge to an edge of T{p); 
W is constructed by traversing T{p) around, with respect to its planar embedding of G. All the 
edges used W are not removed when T{p) is removed, and could not have been removed when some 
previous T{p') was removed. Hence, walk W remains in the graph after T{p) is removed. Since the 
removed dual edges of T{p) always connect two faces visited by W, we infer that all the vertices of 
G* that could have been possibly disconnected by removal of T{p), are actually still connected by W. 
Therefore, after the second step of the construction of TL the graph is still connected. Since removal 
of a degree-1 vertex does not spoil connectivity, we infer that at the end TL is indeed connected. 

We are left with proving the correspondence between objects of F and faces of TL. Since TL is 
a subgraph of G*, every face of TL consists of a union of a nonempty set of original faces of G*, 
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which moreover were connected in G*. Every face of G* corresponds to vertex of G, which means 
that every face f* of Ti contains a nonempty subset of vertices X{f*), which moreover is connected 
in G. Note now that for every p € X, removal of the dual edges of T{p) merges all the faces of 
G* corresponding to vertices of M(p), which means that vertices of every Voronoi region M(p) are 
contained in one face f* G F{%). It remains to show that no face f* G F('H) can accommodate 
more than one Voronoi region. If this was the case, then, since X{f*) is connected in G, there 
would be an edge uv G E{G), u,v G X{f*), that would connect two vertices from different Voronoi 
regions, and whose dual was removed in the construction of Ti. However, in the construction of Ti, 
in the first step we removed only duals to the edges that connect two vertices of the same Voronoi 
region, and in the second step we removed only vertices of degree I, which cannot merge two faces 
of the graph. □ 


Lemma 4.4 shows that the prediagram Lf is a suitable skeleton of the final diagram. More 
precisely, 7i consists of vertices of degree 3 that are connected using long paths consisting of vertices 
of degree 2. From now on we assume that i > 2. If ^ = I then Ti is empty, and if = 2 then H is a 
cycle. Both these cases are degenerate in the argumentation to follow. 

Given the prediagram T-L = we define the Voronoi diagram % = Tijr as follows. To obtain 
7i from Ti, we iteratively take a vertex of degree 2 in the graph, remove it, and add a new edge 
between its neighbors. Thus, all the paths consisting of vertices of degree 2 in T-L are shortened to 
single edges. Note, that in this manner T-L may become a multigraph: 


• There may appear multiple edges between a pair of vertices u and v, in case there was more 
than one path of vertices of degree 2 connecting u and v in T-L. 

• There may appear a loop at some vertex u, in case there was a path of vertices of degree 2 
traveling from u back to itself. 


In the following, we assume that a loop at vertex u contributes with 2 to the degree of u. The 
following lemma formalizes the properties of the Voronoi diagram TL. 

Lemma 4.5. If i > 2, then the diagram TL is a connected 3-regular multigraph embedded on a 
sphere. Moreover, Ti has exactly 2i — A vertices, — 6 edges, and I faces, and the set of faces ofTi 
is exactly the same as the set of faces ofTi. 


Proof. The diagram TL is constructed from TL by taking every path P in TL that begins and ends in 
vertices u, v of degree 3 but travels only through vertices of degree 2, and replacing it with a single 
edge uv (possibly u = v). The embedding of TL into the sphere is defined as follows: The vertices of 
V{Ti) are embedded exactly in the same manner as in the embedding of Ti. Whenever a path P is 
replaced with an edge uv between the endpoints of P, then tl^ embedding of uv is defined as the 
union of the embeddings of the edges of P. Thus, the faces of TL are exactly the same as faces of TL. 

When constructing Ti, we removed all the vertices of degree 2, while all the vertic es o f degree 3 
keep their degrees. Hence TL is 3-regular. Moreover, since TL was connected (Lemma 4.4), then by 
the construction so is Lf. 

The fact that \F(TL)\ = i follows from Lemma 


4.4 


and the fact that the faces of TL are exactly 
the same as faces of TL. The conclusion that |V(?^)| =21 — 4 and \E(TL)\ = 3i — 6 follows from the 
fact that TL is connected and 3-regular, and a simple application of the Euler’s formula, which is 
valid also for connected multigraphs. □ 


In the following, the vertices of V(Lf) will be also call the branching points of diagram TL. The 
following auxiliary lemma shows that loops of the diagram TL have very simple structure. 
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(a) Voronoi partition M and the prediagram T-i 



(b) Diagram T-i 


Figure 7: Example of a Voronoi partition Mj- and the corresponding prediagram and diagram. 
Objects of F are denoted in black. Note that vertices in the prediagram T-Ljr are faces of the original 
graph G. Branching points of the diagram T-Ljr are depicted in yellow, blue, and red, depending 
whether they correspond to type-1, -2, or -3 singular face according to definitions in Section 4.5 


Lemma 4.6. Suppose H is a connected 2>-regular multigraph embedded on a sphere. Then for every 
loop in H, one of the disks in which the loop separates the sphere is a face of H. 

Proof. Consider a loop e at vertex u in H. Since H is 3-regular and loop e contributes with 2 to 
the degree of u, then u is incident to one more edge e'. Then e' is embedded into one of the disks 
into which e partitions the sphere, and by the connectivity of H we infer that the whole rest of H is 
embedded into the same disk. The second disk contains therefore no point of the embedding of H, 
and hence is a face oi H. □ 

4.5 Important faces 

In the algorithm we consider a large family of long sequences of faces of G, in hope of capturing a 
sequence that forms a short balanced separator of the Voronoi diagram imposed by the solution, 
on which the recursion step can be employed. Since balanced separators of this diagram can have 
width as large as we need to take into consideration all the sequences of at most this length. 

Of course, we could consider every face of G as a possible candidate to be used in a separator, but 
then we would have candidate separators, which is too much for the running time promised 

in Theorem 13.11 

In this section we prove the following: one can enumerate in polynomial time a family of at most 
faces of G, called further important faces, such that the candidates for faces used in separators 
can be safely picked from this family. Thus, we arrive at candidates for separators, instead 

of the original , Formally, this section is devoted to the proof of the following result. 

Theorem 4.7. There exists a family £ C F{G) such that 

• \£\ < d^; ^ 

• for every normal subfamily of objects F TV, it holds that VifHj^) C £. 

Moreover, family £ can be constructed in time 0[d‘^ ■ n^«). 
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Let us remark that Theorem 4.7 is needed only for the purpose of reducing the running time 
if we just allowed all the faces of G to be used to build separators, then we would obtain an 
algorithm solving DISJ OINT NETWORK Coverage in time , in many of our 

it actually holds that n = dP^^\ and hence we would obtain the same 


3.2 


applications in Section 

asymptotic running time without the speed-up due to the results of this section. Therefore, the 
result on enumerating a small family of important faces can be considered as optional, and the 
reader interested only in a time algorithm is advised to skip this section. 

We now move on to the proof of Theorem |4.7t which spans the whole this section. However, 
before we proceed, we need some definitions and auxiliary lemmas. Let us fix some normal family of 
objects T" C 2?, and let M = Mj- be the corresponding Voronoi partition. For two distinct objects 
Pi,P2 £ we say that Voronoi regions M(pi) and M(p2) nieet at face / if df contains both a vertex 
of M(pi) and a vertex of M(p 2 )- 

We say that face / is a type-1 singular face for IF, if for a triple of distinct objects pi,P2,P3 £ 
the Voronoi regions M(pi),M(p2),M(p3) pairwise meet at /. Then we say that triple {pi,P 2 ,P 3 ) 
certifies that / is a type-1 singular face. The following lemma shows that there are generally not so 
many type-1 singular faces. 


Lemma 4.8. Suppose {pi,P 2 ,P 3 } = F is a normal family of three objects. Then there are at 
most 2 faces of G that are type-1 singular faces for F. 


Proof. Let Mi = Mj-(pj) for i = 1,2,3. For the sake of contradiction, suppose that there are 3 
different faces fi, f 2 , f 3 , such that Mi, M 2 , all pairwise meet at each of fi, f 2 , f 3 - Construct a 
planar graph H as follows: 


• For each i = 1, 2,3, introduce a vertex Wi inside face fi, and connect Wi to all the vertices of 

dfi. 

• For each j = 1, 2, 3, contract the whole subgraph G[Mj] to a single vertex mj. Note that this 
is possible since G[Mj] is connected by Lemma 

• Remove all the vertices and edges of the graph apart from {wi: i = 1,2, 3}, {mj : j = 1, 2,3}, 
and the edges between them. 


4.2 


Observe that the obtained graph H is isomorphic to After performing the first step of the 

construction we still had a graph with a sphere embedding, since the new vertices and edges could be 
drawn inside faces ft. Then, we applied only edge contractions, vertex deletions and edge deletions. 
Hence, we have found a as a minor of a planar graph, a contradiction. □ 

Now we define the second type of singular faces. We say that face / is a type -2 singular face for 
F, if there exists a triple of distinct objects pi,P2,P3 £ ^ satisfying the following. One vertex vi of 
df belongs to M(pi) and two vertices V2, V3 of df belong to M(p2)- Moreover, if we take a closed 
walk W va. G constructed by concatenating paths in T{p2) from cen(p2) to V2 and V3, and the edge 
V 2 V 3 , then this closed walk separates loc(pi) from loc(p3) on the sphere. Observe that this walk is 
not necessarily a path, but, as it is formed by two paths in the tree T{p 2 ) closed by an edge, its 
removal partitions the sphere into two disks. The requirement is exactly that one of these disks 
contains loc(pi) and the second contains loc(p3). 

Similarly as before, we will say that triple {pi,P 2 ,P 3 ) certifies that / is a type-2 singular face. 
The following lemma shows that, again, there are not so many type-2 singular faces. 

Lemma 4.9. Suppose {pi,P 2 ,P 3 } = F TV is a normal family of three objects. Then there is at 
most 1 face of G that is a type-2 singular face for F and is certified by the triple {pi,P 2 ,P 3 )- 
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(a) Type 1 


Figure 8 : Singular faces of types 
depicted in gray. 


^2 




1, 2, and 3. Paths from a vertex to the nearest object have been 


Proof. Suppose there are two such faces / and f'. We adopt the notation from the definition of a 
type-2 singular face for /: {ui} = V{df) nM(pi), {^2,^3} = V{df) nM(p2), Pt is the cen{p2)-vt 
path inside T{p2), for t = 2,3, and W is the closed walk formed by concatenation of P2, P3, and 
edge V2V3. Note that edge V2V3 does not belong to T{p2), since then P2 would be a subpath of 
P3 or vice versa, and W would degenerate to a path that would not separate any two nonempty 
disks. Let Di be the open disk of S \ kF that contains pi, and D3 be the one that contains P 3 . 
Since vi G M(pi), loc(pi) is contained in Di, and G[M(pi)] is connected, we infer that vi G Di. 
Moreover, if Pi is the cen(pi)-ui path inside T{pi), then Pi is entirely contained in Di. 

Let for t = 1,2, 3, P/ for t = 2,3 and W' be the analogical objects defined for /'. Since /' 
contains a vertex of M(pi), then we have that f is contained in Di. Since P 2 and P 3 are paths 
inside the tree T{p2), and W is defined also by paths inside the tree T{p2), we infer that both P 2 
and P 3 are entirely contained in Di U W. Moreover, they do not contain edge V 2 V 3 (since this edge 
does not belong to T(p 2 )). As /' 7 ^ / and the other face incident to V 2 V 3 lies outside Di, we have 
also that V 2 V 3 7 ^ '^ 2 ^ 3 . Consequently, we infer that the closed walk W' does not use the edge V 2 V 3 . 

To obtain a contradiction, we now exhibit a curve 7 on the sphere that connects cen(pi) with 
cen(p 3 ) and does not touch W'] this will contradict the fact that W’ separates loc(pi) from loc(p 3 ). 
Construct 7 by: 

(i) starting by traveling from cen(pi) to ui using path Pi inside r(pi); 

(ii) traveling inside the face / from vi to the middle of edge V 2 V 3 ; and 

(hi) finishing by traveling inside the disk D3 from the middle of edge V2V3 to cen(p 3 ). 


In part (i) of 7 we do not touch W', because we only travel through vertices of M(pi), whereas all 
the vertices on W' belong to M(p 2 )- Then, in part (ii) we do not touch W, since W does not use 
the edge V 2 V 3 . Finally, in part (iii) we do not touch W' since W is disjoint with D 3 . Hence, 7 is a 
curve from cen(pi) to cen(p 3 ) that avoids W, a contradiction. □ 


Finally, we define the third type of singular faces. We say that face / is a type-3 singular face 
for P, if there exists a quadruple of distinct objects Po,Pi,P 2 ,P 3 £ P satisfying the following. Let 
{ui, U 2 , U 3 } = V(df). The first requirement is that vi,V 2 ,V 3 G M(po); but no edge of df is contained 
in T{pq). Let Pt be the path inside T{po) from cen(po) to vt, for t = 1,2, 3. Let Wt be the closed 
walk formed by concatenation of paths Pt+i, Pi +2 and the edge vt+iVt+ 2 , where the indices behave 
cyclically. Observe that removal of walks Wi, W 2 , W 3 partitions sphere S into four open disks: disks 
Dt for t = 1 , 2 ,3 such that Wt contains the boundary of Dt, and the last disk being simply the face 
/. Then the second requirement is that loc{pt) is entirely contained in Dt, for t = 1,2,3. 


31 





Similarly as before, we will say that a quadruple {po-,Pi,P 2 -,P‘i) certifies that / is a type-3 singular 
face. The following lemma shows that, again, there are not so many type-3 singular faces. 

Lemma 4.10. Suppose {po,Pi,P 2 ,P 3 } = F (fiV is a normal family of four objects. Then there is at 
most 1 face of G that is a type-3 singular face for T and is certified by the quadruple {po,Pi,P 2 ,P 3 )- 

Proof. Suppose there are two such faces / and f'. We adopt the notation from the definition of a 
type-3 singular face for /, and we will follow the same notation but with primes for the face f. By 
assumption, no edge of / or f belongs to T{po). 

Recall that removal of walks Wi, W 2 , W 3 partitions S into four open disks: Di,D 2 , and /. 
Since f f, f is contained in one of the disks Di, D 2 , D 3 . Without loss of generality suppose f is 
contained in Di. Removal of Wi from S partitions S into two disks: one of them is simply Di, and 
the second is a disk Di that contains D 2 , D 3 and /. Observe that since Wi is formed by edge V 2 V 3 
and two paths P 2 ,P 3 , which are contained in the tree T{pq), then each of the paths PfiP!^, P 3 has 
to be entirely contained either in Di U Wi or in Di U Wi — this is because paths PfiP!^, P 3 are also 
contained in the tree T(pq). Since dfi C Di U Wi, we infer that all these paths are contained in 
Di U Wi. Since f is also contained in Di, we conclude that all the closed walks Wfi W'^., W 3 are 
contained in Di U Wi. ^ 

Now examine disk Di. Notice that this disk is disjoint with all the walks W 2 , and does 
not contain f, which means that it is entirely contained in one of the disk D[, D 2 , or D'^. However, 
Di contains both loc(p 2 ) and loc(p 3 ), while every disk D[, D 2 , D'^ contains only one of them. This 
is a contradiction. □ 


In Lemmas 4.8 4.9, and 4.10 we focused only on families J- of size 3, 3, and 4, respectively. The 
following lemma, which follows immediately from Lemma 4.3 and the definition of singular faces of 
types 1, 2, and 3, justifies why we can do it. 


Lemma 4.11. Suppose J-' TP is a normal family of objects. Then for a face f of G, the following 
conditions hold: 


• If f is a type-1 singular face for if, certified by a triple {pi,P 2 ,P 3 ), then it is a type-1 singular 
face for {pi,P 2 ,P 3 }, certified by the same triple ipi,P 2 ,P 3 )- 

• If f is a type-2 singular face for IF, certified by a triple {pi,P 2 ,P 3 ), then it is a type-2 singular 
face for {pi,P2,P3}, certified by the same triple {pi,P2,P3)- 

• If f is a type -3 singular face for IF, certified by a quadruple {po,Pi,P2,P3), then it is a type -3 
singular face for {pq,pi,P 2,P3}, certified by the same quadruple {po,Pi,P2-,P3)- 

Now we show that for a normal family IF, every branching point of FLjr is always a singular face, 
and moreover its type depends on the number of bridges it is incident to in Hjr. Note that since 
lijr is 3-regular, every its vertex is incident to 0, 1, or 3 bridges of Fij^. 

Lemma 4.12. Suppose IF TP is a normal subfamily of objects, and let FL = FLjr be the corresponding 
Voronoi diagram. Then for every branching point f G V{FL), the following holds: 

(1) If f is not incident to any bridge in FI, then f is a type-1 singular face for T. 

(2) If f is incident to exactly one bridge in FI, then f is a type-2 singular face for IF. 

(3) If f is incident to exactly three bridges in Fi, then f is a type-3 singular face for IF. 

In particular, every branching point ofFLisa singular face for IF of one of the three types. 
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Proof. Let e^, e^, 63 be the three edges incident to / in (possibly two of these edges are equal, 
in case / is incident to a loop in Then, in the prediagram T-L, f was incident to three 

edges 63 , 62 , 63 , which are the first edges of the paths contracted to 63 , 62 , 63 , respectively (these 
three edges are pairwise different, since is simple). Edges 63 , 62,63 are dual to edges 61 , 62 , 63 , 
respectively, which form the boundary of the face /. 

We perform a case study depending on how many of edges 63 , 62 , 63 are bridges in Ti (note that 
a loop is not considered a bridge). Observe that, for f = 1 , 2 , 3, edge is a bridge in PL if and only 
if 63 is a bridge in PL. 

Suppose first that none of edges 63 , 62 , 63 are bridges. Then / is incid ent t o three differ ent f aces of 
PL. Since faces of PL correspond one-to-one to objects of J- (see Lemmas 4.4 and Lemmas 4.5), there 


are three distinct objects Pi,P 2 ,P 3 £ -T such that / is incident to corresponding faces f*_^, f*^, f*^ of 
PL. This means that one of the vertices of / belongs to M(p 3 ), second to M(p 2 ), and third to M(p 3 ), 
and the triple {pi,P 2 ,P 3 ) certifies that / is a type -1 singular face for T. This proves [(T)| 

Suppose second that 63 is a bridge in PL, while edges 62,63 are not bridges. This means that 
there exists two distinct objects pi,P 2 £ with corresponding faces of PL being f*^ and f*^, such 
that the face incident to / between 62 , ^3 is /pj, while the face f*^ resides on both sides of the edge 
63 . Let Hi be the subgraph of PL induced by all the branching points that become disconnected 
from / by removing the bridge 63 . Since PL is 3-regular, Hi is not a tree, so it has at least two 
faces. Only one face of Hi contains face f*^, and all the other faces of Hi are actually faces of PL. 
Let then p 3 be any object of LF whose corresponding face f*^ oi PL is a face of Hi. Similarly as 
in the previous case, one vertex vi G V{df) belongs to M(pi) and two vertices of V 2 ,V 3 G V{df) 
belong to M(p 2 )- Moreover these two vertices V 2 ,V 3 are connected by the primal edge ei. Consider 
now a closed walk W obtained by concatenating: walk P 2 inside T{p 2 ) from cen(p 2 ) to V 2 , walk P 3 
inside T{p 2 ) from cen(p 2 ) to V 3 , and edge 61 = V 2 V 3 . Removal of walk W from the sphere separates 
the sphere into two open disks, one of which contains all the branching points of V{Hi) — and in 
particular also loc(p 3 ) — and the second contains all other the branching points of Lf — and in 
particular also loc(p 3 ). Hence, W separates loc(p 3 ) from loc(p 3 ). This means that triple {pi,P 2 ,P 3 ) 


certihes that / is a type-2 singular face for F. This proves (2) 


Finally, suppose that all the edges 63 , 62,63 are bridges in PL. This means that there is one 
object pq £ F and corresponding face f*^ of PL, such that is on both sides of each of the edges 
63 , 62 , 63 . Let {vi,V 2 , U 3 } = V(df), where et = vt+iVt +2 for t = 1,2, 3. Hence, in particular we have 
that vi,V2,V3 G M(po)- Similarly as in the previous case, let Ht be the subgraph of PL induced by 
all the branching points that become disconnected from / by removing the bridge , for t = 1, 2, 3. 
Each Ht is not a forest, since PL is 3-regular. Moreover, each Ht has one face containing the original 
face f*^, whereas all the other faces of PL are partitioned into face sets F{Hi), F{H 2 ), F{H 3 ). 
Hence, there exist three distinct objects Pi,P 2 ,P 3 £ distinct from po, such that the corresponding 
faces /p^, /* 2 , /p 3 of PL are actually faces of Hi, H 2 , and H 3 , respectively. We now dehne Pt to be 
the path inside T{pq) from cen(po) to vt, for t = 1,2, 3, and Wt to be the closed walk formed by 
concatenation of paths Pt+i, Pt +2 and the edge vt+iVt+ 2 - Similarly as in the previous point, we 
see that each walk Wt separates loc{pt) from loc(pf+ 3 ) and \oc{pt+ 2 ). This means that quadruple 


iPo,Pi:P 2 ,P 3 ) certifies that / is a type-3 singular face for F. This proves (3), and concludes the 
proof. □ 


We are finally ready to prove Theorem 4.7, that is, present an enumeration algorithm for 
important faces. Family £ will consist of a union of three families £i,£ 2 ,£ 3 . Families £i and £2 
are constructed by considering every triple objects pi,P 2 ,P 3 £ s.t. {pi,P 2 ,P 3 } is normal, and 
including all the faces that are type-1, resp. type-2, singular faces for {pi,P2,P3}, certified by 
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{pi,P2,P3)- Lemmas 4.8 and 4.9 guarantees that for every triple we have at most two such faces for 


type 1 and at most one such face for type 2, so in total we have that \£i\ < 2d{d — l){d — 2) and 
\£2\ < d{d — l)(d — 2 ). Similarly, family <f3 is constructed by considering every quadruple of objects 
po,pi,P2,P3 s.t. {po,pi,P2,P3} is normal, and including all the faces that are type -3 singular faces 


for {po,pi,P2,P3}, certified by ipo,Pi,P2,P3)- Lemma 4.10 guarantees that for every quadruple we 


have at most one such face, and thus l^lsl < d{d — l){d — 2){d — 3 ). Note that, for a given triple or 
quadruple of objects, it can be verified in polynomial time using the definition whether a given face 
is a type-l/type- 2 /type -3 singular face for this triple/quadruple. Hence, families £1,82, £3 can be 
constructed in time 0{d^ ■ 

We now set <5 = <Si U <?2 U <^3, so in particular we have the following: 

1^1 < 3did - l){d - 2) + did - l)(d - 2){d -3) = d^id- l)(d - 2) < d^. 


Let us consider any normal subfamily of objects T By Lemma [ 4. 11 we have that if / G L’(G) 
is a type-t singular face for T , then / has been included in respective family £t^ for t = 1 , 2 , 3 . On 


the other hand. Lemma 4.12 implies that every branching point of Tijr is a singular face for of 


one of the three types. We infer that every branching point of Tijr belongs to £. This concludes the 


proof of Theorem 4.7 


4.6 Curves, nooses, and sphere cut decompositions 

Curves. Recall that for us a curve on a sphere is a homeomorphic image of an interval, and a 
closed curve is a homeomorphic image of a circle. A directed (closed) curve 7 is simply a (closed) 
curve with a specified direction of traversal. The same curve but with the opposite direction of 
traversal will be denoted by 7“^. For a directed closed curve 7, the area enclosed by 7, denoted 
enc(7), is the open disk of S \ 7 that lies to the right of 7 in the direction of its traversal (in this 
definition we fix the orientation of S). The second open disk, called the area excluded by 7, will be 
denoted by 0x0(7). Note that enc(7) = 0x0(7“^). 


Noos 0 s and th 0 radial graph. Let us consider a situation where some connected multigraph 
H is embedded on the sphere S. We say that a directed closed curve 7 is a noos^ if it intersects 
the embedding of G only in the vertices of H, and moreover visits every face of H at most once. 
The length of a noose 7 is equal to \ ViH) Hql, the number of vertices it traverses; we require that a 
noose traverses at least one vertex. We often say that a noose 7 is with respect to H if the graph H 
is not clear from the context. 

When considering nooses, we shall work on the radial graph of H, denoted Rad(L7). The graph 
Rad(L7) is a bipartite multigraph, where one partite set consists of F(iL), the vertices of H, and 
the second partite set is F(iL), the faces of H. For every face / G T(R) and vertex v G V{H) we 
put an edge between / and v for every occurrence of v on the boundary of /; note that if v appears 
multiple times on the boundary of /, it will be connected to / in Rad(R) via multiple edges. A 
sphere embedding of H naturally induces a sphere embedding of H as follows. For every face / we 
put the corresponding vertex of Rad(R) inside /; we call this vertex the center of /, and denote it 
by C0n(/). Then we connect it to the vertices on the boundary of / radially, according to their 
order of traversal on df. See Figure for an example. 

Now observe that every noose 7 can be homeomorphically transformed to an equivalent noose 
(in terms of the order of vertices and faces visited) that is a simple cycle in the radial graph Rad(L7) 

^In the classic literature nooses are not directed, but it will be convenient for us to work with directed ones. Also, 
often nooses are not required to visit every face at most once and nooses with this property are called tight. For us, all 
the considered nooses are tight. 
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with a chosen orientation. More precisely, whenever 7 travels from a vertex rt to a vertex v through 
a face /, then we replace this part of 7 by going from u to cen(/) along an edge of Rad(i7), and 
then from cen(/) to v again along an edge of Rad(R). Note that if v appears on the boundary 
of / multiple times, we may always pick the edge between cen(/) and v that corresponds to the 
direction from which v is accessed by 7. The same holds for entering / from u. Hence, from now on 
whenever considering some embedded connected multigraph H, we will consider only nooses that 
are oriented cycles in Rad(R). 

Given some noose 7, we may define the subgraph enclosed by 7 as the subgraph enc(7, H) oi 
H consisting of all the vertices and edges of H that are embedded into enc(7) U 7. Similarly, the 
subgraph excluded by 7 is the subgraph 0x0(7, H) oi H that consists of all the vertices and edges of 
H that are embedded into 0x0(7) U 7. Note that (£'(0110(7, R)), £^(0x0(7, R))) is a partition of 
E{H), and 1^(0110(7, £)) n ^(0x0(7, £)) consists of exactly the vertices traversed by 7. 



Figure 9 : The Voronoi diagram from Figure [^together with its radial graph (dotted edges) and an 
exemplary noose 7. The area enclosed by 7 is depicted in green. 


Sphere cut decompositions. We now recall the main tool that will be used to justify the 
correctness of our algorithm, namely sphere cut decompositions. Intuitively speaking, every graph H 
embedded on a sphere has a hierarchical decomposition (formally, a branch decomposition), where 
every separator used to decompose the graph is formed by vertices traversed by some noose of length 
at most 0{^^\V{H)\). In our case we will apply the sphere cut dec omposition to the diagram l-iz, 
where Z is the optimum solution to the problem. By Lemma 4.5 we have that \Vi^z)\ = 0 {k), 


which means that all the separator s use d in the (unknown) sphere cut decomposition of Hz have 
length 0{y/k). By dint of Theorem 4.7 we are able to construct a family of candidates for 

these separators. Then we run a recursive algorithm that at each step guesses a separator that 
splits the (unknown) solution in a balanced way, and then recursively solves simpler subinstances. 
Therefore, to see why the algorithm is correct, it is necessary to have a good understanding of noose 
separators in planar graphs, which brings us to the framework of sphere cut decompositions. 

We now move to a formal introduction of sphere cut decompositions. Recall that for a multigraph 
G, a pair (T, rj) is a branch decomposition of G if T is a tree with all the internal vertices of degree 3 , 
and r/ is a bijection from E{G) to the set of leaves of T. For an edge e G E{T), we define the middle 


35 






















set of e, denoted mid(e), as follows. Suppose that removal of e breaks T into two trees Ti,T 2 , and 
let Fi,F 2 be the sets of edges of G mapped to the leaves of T contained in Ti and T 2 , respectively. 
Then (Fi, F 2 ) is a partition of E{G). Set mid(e) consists of all the vertices of G that are incident 
both to an edge of Fi and to an edge of F 2 . The width of edge e is then defined as |mid(e)|, and the 
width of decomposition (T, rj) is the minimum width among the edges of T. The branchwidth of G, 
denoted bw(G), is the minimum possible width of a branch decomposition of G. 

Suppose now that G is a connected multigraph without loops that is embedded on a sphere S. 
A sphere cut decomposition {sc-decomposition) of G is a triple {T,r],5), where: 


• (r, r]) is a branch decomposition of G; 

• (5 is a function that assigns to every edge e G E{T) a noose (5(e) that traverses exactly the 
set of vertices mid(e). Moreover, if (Fi,F 2 ) is the partition of E{G) as in the definition of 
mid(e), and Gi, G 2 are the subgraphs of G spanned by Fi, F 2 , respectively, then one of Gi, G 2 
is enc((5(e),G) and the second is exc(5(e),G). 

Note that in the second condition, whether enc((5(e),G) is Gi or G 2 depends on the choice of 
direction of the noose (5(e). Observe also that the given partition (Fi,F 2 ) of F(G) uniquely defines 
the noose (5(e) (treated as a simple cycle in Rad(G), without orientation), for the following reason: 
For every face / that is traversed by (5(e), the intersections of df with Fi and F 2 have to be just 
two intervals of the walk df, since face / is visited only once by <5(e). Hence, by examining these 
intervals we can uniquely deduce between which vertices the noose travels within /, and from which 
directions these vertices are accessed. 

We will regard sc-decompositions also as branch decompositions, with all the inherited notions 
for them. 

As observed by Dorn et al. |19j . the following theorem was implicitly proved by Seymour and 
Thomas in [l5]. The improved running time (from original O(n^) to 0{n^)) is due to Gu and 
Tamaki [26]. 


Theorem 4.13 t [im [261145] 1. Suppose G is a bridgeless, connected multigraph G without loops, 
embedded on a sphere S. Then there exists an sc-decomposition {T,r],5) such that {T,r]) is a branch 
decomposition of G of optimum width. Moreover, such an sc-decomposition can be found in time 
G(ra3). 


Let us now deliberate on the proof of Theorem 4.13 as the statement differs in important details 
from the one provided by Dorn et al. |19) (see Theorem 1 in this work). The following argumentation 
is based on the sketch of the proof of Theorem 1 in m- 

The proof follows from (5.1) and the proof of (7.2) in |l5]. In (5.1), it is proved that for every 
2-connected multigraph G, there exists an optimum carving decomposition of G that uses only bond 
carvings (is a bond carving decomposition). Here, a carving decomposition is a similar concept to a 
branch decomposition, but this time leaves of decomposition T correspond one-to-one to vertices of 
G, and the width of an edge e G E{T) is defined as the number of edges traversing the corresponding 
partition (Ai, A 2 ) of V{G). This partition (Ai, A 2 ) is called a carving, and the carving is a bond 
carving if both Ai and A 2 induce connected graphs in G. Equivalently, the edges traversing the 
carving form an inclusion-wise minimal edge cut between some two vertices of G. 

Note now that if G is embedded on a sphere S, then such a minimal cut corresponds to a simple 
cycle in the dual G* , which traverses duals of the consecutive edges of the cut. This cycle divides 
the sphere into two open disks such that vertices of Ai are embedded into one of them and vertices 
of A 2 are embedded into the second of them. 

Carving decompositions are then translated to branch decompositions using the notion of a 
medial graph Med(G), which is just the dual of the radial graph. More precisely, vertices of the 
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Med(G) are the edges of the primal graph, and for every vertex v £ V(G) we connect the edges 
incident to v into a cycle, in the order of their appearance around v in the sphere embedding. Thus, 
the medial graph is also embedded into the same sphere. In the proof of (7.2) in |45j . Seymour 
and Thomas consider carving decompositions of the medial graph Med(G). Since the vertex set 
of Med(G) is exactly the edge set of G, carving decompositions of Med(G) correspond exactly 
to branch decompositions of G. Moreover, cycles in the dual of Med(G) (which is just Rad(G)) 
correspond exactly to nooses in G. Hence, bond carving decompositions of Med(G') correspond 
to sc-decompositions of G. Seymour and Thomas prove that this correspondence translates an 
optimum bond carving decomposition of Med(G) into an optimum branch decomposition of G, 
where the width gets divided by two. In particular the carving width of Med(G) is twice larger 
than the branchwidth of G. Since this translation produces in fact an sc-decomposition, we conclude 
that there exists an sc-decomposition of G of optimum width. 

In (9.1) Seymour and Thomas give an 0{n‘^) algorithm that constructs an optimum bond carving 
decomposition of an planar graph. By applying this algorithm to Med(G) and translating the 
output to an optimum sc-decomposition of G, we obtain the sought result. The running time of the 
algorithm has been improved from original 0{n‘^) to O(n^) by Gu and Tamaki |26j . 

Observe now that in (5.1) of [IS], it was assumed that the graph, whose carving decompositions are 
considered, needs to be 2-connected. Indeed, for every vertex v £V {G) we have that ({u}, V (G)\{u}) 
is one of the carvings in every carving decomposition. If now v was a cutvertex, then this carving 
would not be a bond carving. 

It is now easy to see that the medial graph Med(G) is 2-connected if and only if the primal 
graph G is connected and bridgeless. Therefore, the whole presented argumentation holds provided 
that G is connected and bridgeless. And indeed, if G contained some bridge b, then {{b},E{G) \ {6}) 
would be one of the partitions {Fi,F 2 ) considered for some edge of the decomposition T. However, 
then the noose enclosing only the bridge b would need to visit the face on both sides of b twice. 


which contradicts the definition of a noose. Therefore, in Theorem 4.13 it is necessary and sufficient 
to assume that G is connected and bridgeless. 

For the other assumptions, Dorn et al. m consider only simple graphs. However, all the 
arguments work in the same way in presence of multiple edges, since the original proof of Seymour 
and Thomas |l5] works for multigraphs (as is noted in the first sentence of [35]). We need, however, 
to exclude loops, since in the presence of loops, understood as curves connecting a vertex to itself, 
there arise technical problems with the medial graph Med(G). Namely, for the natural generalization 
of this notion to multigraphs with loops, a loop can correspond to a cutvertex in Med(G). We 
remark that Seymour and Thomas allow loops in their work, but they understand them differently, 
as edges of arity I, being de facto annotations at vertices. We also remark that in their version of 


Theorem 4.13 Dorn et al. m (incorrectly) assume only that the graph does not have vertices of 
degree 1, instead of excluding bridges. As we have seen, this assumption is necessary. 

The second necessary ingredient is the well-known fact that planar graphs have branchwidth 
bounded by approximately the square root of the number of vertices. 

Theorem 4.14 (|23j). For a planar graph G on n vertices, it holds that bw(G) < \/F5n. 


We note that Theorem 4.14 can be trivially extended to multigraphs. Thus, Theorems 4.13 


and 4.14 together imply that every bridgeless, connected multigraph G without loops that is 
embedded on a sphere, admits an sc-decomposition of width at most y^4.5|I7(G)|. This corollary 
will be our main tool in the sequel. 
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4.7 Balanced nooses in plane graphs 

Given a connected multigraph G embedded on a sphere, a noose 7 w.r.t. G is a-edge-balanced if the 
numbers of edges enclosed and excluded by 7 are bounded by a\E{G)\, i.e. |i?(enc( 7 , G))| < a\E{G)\ 
and |-B(exc( 7 , G))| < a\E(G)\. The following lemma shows that the plane multigraphs we are 
interested in admit short edge-balanced nooses. 


Theorem 4.15. Let G be a connected 3-regular multigraph with n vertices, m > 6 edges, possibly 
with loops, and embedded on a sphere S. Then there exists a '^-edge-balanced noose w.r.t. G that 
has length at most \/A.5n. 

Proof. Let T* be the tree of bridgeless components of G, where with every node x G V(T*) we 
associate a subset C V{G) inducing a bridgeless component of G, and every edge xy G E{T*) 
corresponds to a bridge bxy connecting G[Bx] with G[By]. For every edge xy G E{T*), consider 
the connected components Gi and G 2 of G — bxy and let us direct this edge from x to y if 
|ili(Gi)| > \E{G 2 )\, from y to x if |i?(Gi)| < |£'(G 2 )|, and breaking the tie arbitrarily in case 
|Fi(Gi)| = \E{G 2 )\. Thus, every edge of T* becomes directed. Since T* has exactly |F(T*)| — 1 
edges, we have that the sum of outdegrees in T* is equal to |F(T*)| — 1, and hence there exists a 
node xo G V{T*) that has no outgoing directed edge. 

Let H = G[Bxq] and observe that H is a, bridgeless, connected multigraph embedded on S. 
Moreover, observe that H does not contain loops for the following reason: A loop in a 3-regular 
multigraph must be attached to a vertex having exactly one other incident edge that is a bridge. 
Hence, if H contained a loop then Bxq = V{H) would consist of one vertex and xq would be a leaf 
in T*. Since the only edge incident to xq in T* is directed towards xq and G[Bxq] has one edge, we 
would have that m<2-l-(-l = 3. This would be a contradiction with the assumption that m > 6 . 

Suppose first that \Bxq\ = 1, i.e., H is graph consisting of a single vertex u. Then u is incident to 
3 bridges in G, and removal of each of these bridges separates from u a connected component having 
at most edges. Let G be the component that has the most vertices among these three ones; 
then \E{C)\ > = y — 1. Take a noose 7 of length 1 that visits u and the unique face incident 

to u, chosen in such a manner that it encloses exactly component G and the bridge connecting G 
with u. Then, as m > 6 , it follows that 


|.F(enc( 7 ,G))| 

|F;(exc( 7 ,G))| 


, „ / m — 1 2 

\E{G)\ + 1 < + 1 < -m, 

/ TTi \ 2 

m - {\E{G)\ + 1) <m - - 1 1 j = -m. 


Hence, 7 satisfies the required properties. 

From now on suppose that \Bxq\ > 1. We apply Theorems 4.13 and 4.14 to H in order to find 
an sc-decomposition (T, y, 5) of H of width at most y^4.5\V{H)\ < \/4.5n. Define a weight function 
w: E(H) —>■ N as follows: We first put w(e) = 1 for all e G E[H). Then, for every bridge uv such 
that u G Bxq and v ^ Bxq, consider the connected component G oi G — uv that does not contain 
Bxg and add its edge count (including the bridge uv itself) to the weight of one of the edges of H 
incident to u; such an edge always exists for H is bridgeless and has more than one vertex. Moreover, 
every vertex of H is incident either to three edges of H or to two edges of H and one bridge, and 
hence it is easy to see that the distribution of weights can be done in such a manner that to the 
weight of every edge e of Lf we add the weight corresponding to at most one bridge incident to an 
endpoint of e. Observe that thus we have that w{H) = m, where the weight of a graph is defined 
as the sum of the weights of its edges. Moreover, since G is 3-regular and every edge of T* incident 
to Xq was directed towards xq, we have that w(e) < 2 -|- = y | for each e G E{H). 
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Now, consider any edge e G E(T). Let {Fi, F 2 ) is the partition of E{H) as in the definition of 
mid(e), and let Lfi, H 2 be the subgraphs of H spanned by Fi, F 2 ; hence we have that one of Ffi, H 2 
is enc{5{e), H) and the second is exc(h(e), iL), depending on the orientation of <5(e). For every 
e G E{T), let us select this orientation so that w(enc((5(e), i^)) < w(exc((5(e), FT)), breaking ties 
arbitrarily in case of equality. Also, let us direct the edge e in T so that it points to the subtree of T 
whose leaves span exc(S(e), G). Again, since after directing edges of T the total sum of outdegrees 
in |l^(r)| — 1, there exists a node 2 ; G y{T) such that there is no outgoing edge from z in T. 

Consider first the case when 2 ; is a leaf in T; let uv = rj{z) and let ez be the unique edge of 
T incident to z. Hence, 5{ez) is a noose that excludes only the edge e^, and encloses the whole 
rest of E{H). Since was directed towards z, we have that w{uv) = w(exc(h(e^),FF)) > y. 
Since m > 6, this means that one of the endpoints of uv, say u, is incident to a bridge uw such 
that the weight corresponding to this bridge contributed to the weight of uv. More precisely, if 
C is the connected component of G — uw that does not contain Bxq, then ■w(ttu) was increased 
from initial value of 1 by contribution 1 + \E{C)\. Let 7 be the length-1 noose that visits u and 
the face around bridge uw, and encloses exactly E{C) and the bridge uw. Then we have that 
|FF(enc( 7 , G))| = \E{C)\ + 1 = w{uv) — l<y + § — 1< |m, as m > 6. On the other hand, 
|FF(exc( 7 , G))| = m — (|£'(C')| -|- 1) = m — (w(uu) — l)<m—y-|-l< again as m > 6. Hence, 
7 satisfies the required properties. 

Finally, we are left with the case when 2; is an internal vertex of T. Let ei, 62,63 be the edges of 
T incident to 2;; recall that they are all directed towards 2;. For t = 1 , 2,3, let Ht = enc(<5(et), FF). 
Then we have that w(FFi) + w(FF 2 ) + w(FF 3 ) = w(FF) = m, and < y for all t = 1,2,3. 

W.l.o.g. suppose that w(FFi) is the largest among w(FFt) for t = 1,2,3; then w(FFi) > y. Let 
7 ' = S{ei), which is a noose w.r.t. FF; recall that the length of ^ is at most y^ 4.5|1F(FF)| < \/4.5n. 
Let us create a noose 7 w.r.t. G by modifying noose 7 ' as follows: whenever traverses a vertex 
u € V (FF) that in G is incident to a bridge uw connecting u with a component C, then draw 7 
around u in such a manner that it leaves uw together with the whole component G on the same 
side of 7 as the edge of FF whose weight G was contributing to (see Figure 10). Note that since 
G is 3-regular, this is always possible. In this manner, we have that |FF(enc( 7 , G))| = w(FFi) and 
|FF(exc( 7 , G))| = w(FF 2 ) -|- w(FF 3 ), and 7 has the same length as 7 '. 



Figure 10: The modification applied to construct 7 from 7 ' around vertices of Bx^ that are incident 
to a bridge in G. The left panel shows noose w.r.t. FF, and the right panel shows the adjusted 
noose 7 after reintroducing component G. 


Of course, since w(FFi) < y, we have that |FF(enc( 7 , G))| = w(FFi) < y < |m. On the other 
hand, since w(FFi) > y, we have that |Fi'(exc( 7 , G))| = m — w(FFi) < |m. Thus, 7 satisfies all the 
required properties. □ 

For the sake of later argumentation, we need a slightly different definition of balanceness that 
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concerns faces instead of edges. We say that a noose 7 is a-face-balanced if the numbers of faces of 
G that are strictly enclosed and strictly excluded by 7 , respectively, are not larger than a|F(G)|. 
Here, F(G) is the set of faces of G. Note that faces traversed by 7 are neither strictly enclosed nor 
strictly excluded by 7 , and hence they contribute to neither of these numbers. 

Lemma 4.16. Let G be a connected 3-regular multigraph, possibly with loops, and embedded on a 
sphere E. If a noose 7 w.r.t. G is “^-edge-balanced, then it is also ^-face-balanced. 

Proof. Since G is connected and 3-regular, by Euler’s formula we have that |E(G)| —|iil(G)|-|-|F(G)| = 
2 and by degree count we have that |E^(G)| = ||E(G)|. From this it follows that 

\F{G)\=2 + ^\E{G)\. (1) 

As G is connected and noose 7 visits every face of G at most once, it follows that both 
Gi := enc(G, 7 ) and G 2 ■= exc(G, 7 ) are connected. Hence, Gi is a connected plane multigraph, 
possibly with loops, with maximum degree at most 3. From Euler’s formula it follows that 
|l/(Gi)| — |Fi(Gi)| -|- |F(Gi)| = 2 , and from the degree bound we infer that |F(Gi)| < ||F(Gi)|. 
Therefore, 

|F(Gi)| = 2-|y(Gi)| + |F(Gi)|<2-?|F(Gi)| + |F(Gi)| = 2 + l|F(Gi)| 


here, the second inequality follows from the assumption that 7 is |-edge-balanced, whereas the last 
equality follows from Q. Observe now that the number of faces of G strictly enclosed by 7 is equal 
to |F(Gi)| — 1; indeed the face set of Gi comprises exactly the faces of G strictly enclosed by 7 plus 
one extra new face that contains the curve 7 . Hence, the number of faces of G strictly enclosed by 
7 is at most | -|- ||F(G)| — 1 < ||F(G)|. A symmetric argument shows that also the number of 
faces strictly excluded by 7 is at most ||F(G)|. □ 


From Theorem 4.15 and Lemma 4.16| we infer the following corollary. 


Corollary 4.17. Let G be a connected 3-regular multigraph with n vertices, m > 6 edges, possibly 
with loops, and embedded on a sphere U. Then there exists a ^-face-balanced noose w.r.t. G that 
has length at most ^/A.5n. 


4.8 Voronoi separators 

Armed with a good understanding of Voronoi diagrams and knowledge about the existence of short 
balanced nooses in plane graphs, we can combine these two ingredients to design a Divide&Conquer 
algorithm for DISJOINT Network Coverage. More precisely, if Z is the optimum solution to the 
considered instance of Disjoint Network Coverage, then the algorithm will iterate through all 
the possible candidates for balanced nooses of the Voronoi diagram TLz ■ Each candidate noose will 
separate the instance into a number of subinstances, which will be solved recursively. For the correct 
selection of the balanced noose of TLz, each of the subinstances will contain only at most objects 
from the solution, and hence we will apply the algorithm recursively only to instances with at most 
this value of the parameter; This will ensure that the running time is as promised in Theorem]^ 
However, first we need to understand formally what are the “possible candidates” for balanced 
nooses, and in what sense they separate the instance at hand into subinstances; these two questions 
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are the topics of the current and the following section. More precisely, we will investigate the 
properties of Voronoi separators, which are structures to be used as separators in the forthcoming 
Divide&Conquer algorithm. 

Whenever we have some normal subfamily of objects J- ^T), then we have a corresponding 
Voronoi partition Mj- and Voronoi diagram T-Lj^. The intuition now is that the nooses w.r.t. Tijr 
can be projected back to the original graph G. That is, suppose in some noose 7 we travel from 
some center of a face cen(/) to an adjacent vertex u. Then in G this corresponds to traveling from 
the center cen(p) of the corresponding object p to the corresponding branching point of T-Ljr using a 
path inside the tree Tjr[p). The crucial observation now is that, for a noose 7 , the knowledge of the 
order in which 7 visits faces of T-Ljr (corresponding to objects of F) and branching points of Fj^, is 
sufficient to deduce the whole projection to G. This is because when projecting we can simply use 
shortest paths. 

Voronoi separators. We move on to a formal argumentation. A Voronoi separator of length r is 
a sequence S of the following form: 


S={ pi,ui,fi,vi, 
P2,U2,f2,V2, 


Pr,Ur,fr,Vr ) (2) 

For this sequence, we put the following requirements (from now on, indices behave cyclically): 

(a) pi,p 2 , ■ ■ ■ ,Pr are pairwise disjoint objects from V, forming a normal family 'D{S) = {pi,p 2 ,...,Pr}- 

(b) /i, / 2 , • • • j /r are pairwise different faces of G. 

(c) For t = 1 , 2 ,... ,r, ut and vt are two different vertices on face ft. Moreover, vt-i,ut G M.'p(^s-^{pt), 
for all t = 1,2,..., r. In other words, vt-i and ut are closest to the location of pt in terms 
of the distance measure (p{v,p) = dist(u, loc(p)) — r{p), among all the locations of objects 
traversed by S. 


We shall treat two Voronoi separators that differ only in a cyclic shift of the sequence as the same 
separator. However, similarly as with nooses, reverting the sequence results in a different separator, 
which we shall denote by S~^. Observe that, given a sequence like in Q, we can verify in polynomial 
time whether it satisfies all the required properties of a Voronoi separator. 

For t = 1,2,... ,r, let Pt be the path from cen{pt) to ut obtained by concatenating a path 
inside T{pt) from cen{pt) to the vertex of loc(pi) closest to ut, with the shortest path from ut to 
loc(pi). Similarly define Qt+i for vt and pt+i. By property (c), paths Pt and Qt are subpaths of 
the tree Txi(^s'^{pt), for every t = 1,2 ,... ,t. Consequently, paths {Pt, Qt}t=i, 2 ,...,r are all pairwise 
vertex-disjoint, apart from paths Pt and Qt that can intersect on a common prefix from the side of 
cen(pt), and otherwise they are vertex-disjoint. 

For a Voronoi separator S, the perimeter oi S is defined as F(S') = IJi=i ^i^t) U V{Qt), i.e., it 
is the union of the vertex sets of all the paths Pt and Qt. Moreover, we define a directed closed 
walk W(S) in G by concatenating consecutive paths Qi, Pi,Q 2 , P 2 , ■ ■ ■, Qt, Pt, where Pt and Qt+i 
are joined using edge utVt, which lies on the boundary of face ft. Note that V{W{S)) = r(5). See 
Figure [TT] for an example. 

Note that, by the definition of a Voronoi separator, VF(5) is a simple cycle in G, possibly with 
simple paths attached to different vertices that correspond to common prefixes of paths Pt and Qt. 
In particular, it still holds that removal of VF(S') from S partitions S into two open disks, one on 
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Figure 11: An exemplary Voronoi separator S of length 4, together with the corresponding walk 
1F(5) (blue). The Voronoi regions of Mx)( 5 ) have been depicted in respective light colors. 


the right side of W{S) (denoted enc(5)), and the second on the left side (denoted exc(S')). For 
any object p G 'D\ 'D{S), we shall say that p is strictly enclosed by S if all the vertices of loc(p) are 
embedded into enc(5) (in particular loc(p) n r(S') = 0). Similarly, p gV\ T){S) is strictly excluded 
by S if all the vertices of loc(p) are embedded into exc(5). Finally, a client q G C is enclosed by S, 
resp. excluded by 5, if pla(g) is embedded into enc(S'), resp. into exc(5); Note that thus every 
client q such that pla(g) G r(S') is both enclosed and excluded by S. Observe also that any path 
in G that connects the location of a strictly enclosed object with the location of another, strictly 
excluded object, must necessarily cross walk kF(5), so it has to traverse a vertex of r(S'). The 
same observation holds also for paths connecting the location of a strictly enclosed object with the 
placement of an excluded client, and vice versa. 


Separators inherited from the diagram. Assume now that we have a normal subfamily of 
objects T" C T>, and let M = Mj- and % = Idj^ be the corresponding Voronoi partition and diagram. 
Suppose further that 7 is a noose w.r.t. T-L. We dehne the corresponding Voronoi separator 5 '( 7 ) as 
follows; see Figure for a visualization. 

Suppose that 7 visits faces and branching points /i, /i, /I, / 2 ) • • • j /r , /r in this order, where 
are different faces of T-L, and ft are different branc hing poi nts o f TL, that is, faces of G (recall that Ti 


is constructed from the dual of G). By Lemmas 4.4 and 4.5, faces of Ti correspond one-to-one to 
objects of F. Let pi,p 2 , • ■ • ,Pr £ T" be the objects that correspond to /i, /Ij • • •, /r j respectively. 
Suppose that when entering ft from ff, for some t G {1,2,... ,r}, noose 7 entered it between edges 
e{ and e^. Let e{ and 62 be the edges of the prediagram TLjr incident to ft that got contracted 
onto el and e^ respectively; in case el or e^ is a loop, we choose the edge of Tijr corresponding to 
the endpoint of the loop by which ft was entered from f^. Then we define ut to be the common 
endpoint of the primal edges ei and 62 corresponding to el and e^, respectively (recall that / is 
a triangle, since G is triangulated). We define vt analogically, based on how 7 leaves branching 
point ft to face ff+i- Note that since a noose, when crossing some branching point, never leaves this 
branching point to the same face and using the same direction as it entered, we have that ut 7 ^ vt for 
alH = 1, 2,..., r. Then separator ^( 7 ) is dehned using pt, ft, ut, and vt as in formula Q. Note that 
£'( 7 ) has the same length as 7 , and moreover, = 5 ( 7 “^). The following lemma encapsulates 

the main properties of separators inherited from the diagram that we shall use later on. 


Lemma 4.18. If j is a noose w.r.t. Ti, then: 
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Figure 12: Separator S{^) and corresponding walk W{S{j)) for the noose 7 from Figure 


(i) S := S{'y) is a Voronoi separator; 

(n) r(5)CUp6C(s)M^(p); 

(Hi) For every object p £ F such that the face ofH corresponding to p is strictly enclosed by 7 
(resp. strictly excluded by ^), we have that p is strictly enclosed by S (resp. strictly excluded 
by S). 


Proof. For (i), the only non-trivial condition is the last one, that Ut,Vt-i G Mx)(s')(pt 
t = 1,2,... ,r. However, from the construction of ut and vt-i and properties of the Voronoi diagram 
(see Lemmas |4.4| and 4.5), it follows that ut,vt-i G 


^APt 


ut,vt-i G M 2 p( 5 )(pi) follows from Lemma 


m 


UpeD(s) • 


). As loc(pt) C 
A{p)- 


for every 
L diagram 

Since 'D{S) C F, the fact that also 
.j^ipt) for all t = 1 , 2 

A 

r{pt), we have that V{Qt),V{Pt) V Mj-(pt), and hence r(S') C 


4.3 


r. 


For (ii), in the paragraph above we have argued that ut,vt-i G 
By Lemma 4.2, also the unique shortest paths from ut and vt-i to loc(pt) are entirely contained 
Apt 


Property (iii) follows directly from (ii) and the construction of the Voronoi diagram T-L. □ 


4.9 Interaction graph and separability of the problem 

Suppose {G,V,C,k) is an instance of Disjoint Network Coverage. By the interaction graph 
C = £(G,2?,C, k) we mean a graph with vertex set VUC, where (a) for p G V and q G C we put 
pq G E{C) if and only if object p covers client q, and (b) for p,p' gV we put pp' G E{C) if and only 
if p and p' contradict the normality requirement, i.e., family {p,p'} is not normal. 

We now define the final abstraction of a separator, which we call a guarded separator. A guarded 
separator is simply a pair X = (Q, F), where Q C P is a normal subfamily of objects and F C V{G) 
is an arbitrary subset of vertices. We will use operators Q(X) and r(X) to extract the first and the 
second coordinate of pair X, respectively. The length of a guarded separator is simply |Q(X)|. Note 
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that every Voronoi separator S naturally induces a guarded separator X(S') = (D(S'),r(S')); this 
will be the main source of guarded separators in our algorithm. 

Let X = (Q, r) be a guarded separator. We say that a client q £ C is covered by X if (7 is covered 
by an object belonging to Q. We also say that an object p £ 'D\Q is banned by X if either Q U {p} 
is not normal, or there exists u G L such that dist(u, loc(p)) — r{p) < dist(u, loc(p')) — r{p'), where 
p' is the object of Q closest to v in terms of the distance measure 4 >{v,p') = dist(u, loc(p')) — r{p'). 
An object p £V\Q that is not banned by X is said to be allowed by p. By cov(X) we denote the 
set of clients covered by X, and by ban(X) the set of objects banned by X. Obviously, by Lemma 4.1 
every object p £V\Q such that loc(p) n L 7 ^ 0 is banned by X. 

We say that two allowed objects p,p' £ V \ Q are separated by X if the shortest path in G 
between loc(p) and loc(p') traverses a vertex of r(X). Similarly, an allowed object p £ 'D\Q and a 
client q £ C are separated by X if the shortest path between loc(p) and pla((7) traverses a vertex of 

r(x). 

The intuition is that if a Voronoi separator S is inherited from some solution, then the objects 
banned by X(S') for sure are not used by the solution. This is because the inclusion of any banned 
object would either contradict the normality, or the fact that walk W (S') traverses only the Voronoi 
regions corresponding to T>{S) in the diagram induced by the whole solution. Thus, the clients 
covered by S and the object banned by S form a “border” that separates the part of the solution 
enclosed by S from the part excluded by S. Here, by “separates” we mean the definition of separation 
for the guarded separator X(S). In the following two lemmas we formalize this separation property. 

Lemma 4 . 19 . Suppose X = (Q, T) is a guarded separator. Suppose further that an allowed object 
p £V\Q and a client q £ C are separated by X. Then the following implication holds: if p covers q, 
then q is covered by X. 

Proof. Suppose p covers q, and let P be the shortest path in G between loc(p) and pla(g). By the 
dehnition of separation, we have that P traverses a vertex of T, say v. Let p' be the object that is 
closest to V among objects of Q in terms of (j){v,p'). Since p is allowed by X, we have that 


dist(u, loc(p)) — r{p) > dist(u, loc(p')) — r(p'). 


( 3 ) 


Since p covers q, we also have that 


dist(pla(g), loc(p)) < s{q) + r{p). 

Finally, since P is the shortest path between loc(p) and pla(q), we have that 


dist(pla(g), loc(p)) = dist(u, loc(p)) + dist(pla(g), u). 


Using ([ 4 ]), and ([^, we infer that 


dist(pla(g), loc(p')) < 
< 


< 


dist(u, pla(g)) + dist(u, loc(p')) 
dist(u, pla(g)) + dist(u, loc(p)) + r{p) 
dist(pla((jr), loc(p)) + r(p') — r{p) 
s{q) + r{p). 


r{p) 


( 4 ) 

( 5 ) 


This means that client q is covered by p' £ Q. □ 

Lemma 4 . 20 . Suppose X = (Q, T) is a guarded separator. Suppose further that allowed objects 
Pi-iP 2 £'D\ Q are separated by X. Then the family {pi,P2} is normal. 
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Proof. Let P be the shortest path in G between loc(pi) and loc(p 2 )- By the definition of separation, 
we have that P traverses some vertex of L, say v. Let p' be the object of Q that is closest to v, in 


terms of (j){v,p'). Since both pi and p 2 are allowed by X, we have that 

dist(loc(pi),n) + dist(loc(p'),?;) > dist(loc(pi),loc(p')) > \r{pi) — r{p)\ (6) 

dist(loc(p 2 )) + dist(loc(p^), n) > dist(loc(p 2 ),loc(p )) > |r(p 2 ) - r{p)\ (7) 

dist(loc(pi), n) — r{pi) > dist(loc(p'), n) — r(p') (8) 

dist(loc(p 2 ),'i') — r{p 2 ) > dist(loc(p'),u) — r(p'). (9) 

Since P is the shortest path between loc(pi) and loc(p 2 ); we have 


dist(loc(pi), loc(p 2 )) = dist(n, loc(pi)) + dist(t', loc(p 2 ))- (10) 

By adding 0 and Q, we obtain: 

dist(loc(pi), ti) + dist(loc(p 2 ), v) + dist(loc(p'), ti) — r(pi) 

> k(P 2 ) — r{p^)\ + dist(loc(p'), v) — r{p') 

> dist(loc(p'), ?;) — r(p2)- 


Thus, by (10) we obtain that 


dist(loc(pi),loc(p2)) > r{pi) -r{p 2 ). 


By adding equations Q and Q, we obtain: 

dist(loc(pi), u) + dist(loc(p 2 ), v) + dist(loc(p'), u) — r{p 2 ) 

> \r{pi) — r{p)\ + dist(loc(p'), v) — r{p') 

> dist(loc(p'), u) — r(pi) 


Again by (10) we obtain that 


dist(loc(pi),loc(p2)) > r{p 2 ) -r{pi). 


( 11 ) 


( 12 ) 


Equations (© and ( [I^ together imply that dist(loc(pi), loc(p 2 )) > V{Pi) — ''’( 7 * 2)15 which means 
that {pi,P 2 } is normal. □ 


For a guarded separator X, we will consider the modified interaction graph >C(X) = C — (Q(X) U 
ban(X) U cov(X)). Let T" C 2? be a normal subfamily of objects, and let M = Mj- and PL = PLjr 
be the corresponding Voronoi partition and diagram. We will say that a guarded separator X is 
compatible with T if Q(X) C LP and r(X) C IJpeQ(x) ®1I.?='(7')' We will moreover say that a guarded 
separator X compatible with P is a-interaction-balanced w.r.t. P if every connected component 
of >C(X) contains at most a\P\ objects from P. Observe that if X is compatible with T", then 
P n ban(X) = 0 by the definitions of banned objects and compatibility. 

Lemma 4.21. If ^ is an a-face-balanced noose in PL, then X(S'( 7 )) is compatible and a-interaction- 
balanced w.r.t. P. 


Proof. Let S = S{j) and X = X(S'). Lemma 4.18 ([(I)] and [(ii)] ) implies that X is compatible with P, 
so by the definition of banned objects we have that P n ban(X) = 0. Hence P\ V{S) C 1/(£(X)). 

Observe that whenever we have two allowed objects pi,p 2 G iF\ V{S) such that pi is strictly 
enclosed by S and p 2 is strictly excluded by S, then loc(pi) and loc(p 2 ) he entirely in two different 


45 






connected components of G — r(5), since any path connecting loc(pi) and loc(p 2 ) has to cross 
the walk Vk(S'). Therefore, any such objects pi,P 2 are separated by X, and Lemma [4.20 implies 
that they are not adjacent in £(X). Similarly, if an allowed object p £ T \ 'D{S) is strictly enclosed 
by S and a client q £ C is excluded by S (or vice versa), then they also lie in different connected 
components of G — r(S') (or q lies on r(S')), and hence are separated by X. Then Lemma 4.19 
implies that p and q are not adjacent in £(X). 

Concluding, every connected component of £(X) either consists only of objects strictly enclosed 
by S a nd clients enclosed by S, or of objects strictly excluded by S and clients excluded by S. 
Lemma 


4.18 


(hi) ensures that objects of T" corresponding to faces of H strictly enclosed (resp. 


excluded) by 7 are exactly the objects of T" that are strictly enclosed (resp. excluded) by S. Since 
7 is a-face-balanced, we infer that there are at most a\F\ objects of T that are strictly enclosed by 
5, and the same holds also for objects strictly excluded by T. Hence, every connected component 
of £(X) can contain at most a\J- \ objects of T", and X is a-interaction-balanced. □ 


We are finally ready to prove the main result of this section, that is, an enumeration algorithm 
for candidates for balanced guarded separators. This is the result whose simplified variant was 
Lemma O 


Theorem 4.22. There exists an algorithm that, given an instance I = {G,'D,C,k) of Disjoint 
Network Coverage with > 4 , enumerates a family Af of guarded separators with the following 
properties: 

(i) |AA| < (24)^®^; and 

(a) for every normal subfamily if TV of cardinality exactly k, there exists a guarded separator 
X G AA that is ‘^-interaction-balanced for T. 

The algorithm works in total time ( 2 d)^®^ • (dcnj^^^^ and outputs the guarded separators of Af one 
by one, using additional (working tape) space {dcn)^^^\ 

Proof. We first apply the algorithm of Theorem |4.7| to compute a family £ of important faces of size 
at most Then, we enumerate all candidates sequences for Voronoi separators of length at most 
3y/k as in formula Q, where faces ft are chosen from the family £. Note that for a fixed length r 
(1 < r < 3Vk), we have at most dJ' ■ d‘^^ ■ 6 ^ such candidates, since every object pt is chosen among 
d options, every face ft is chosen among at most options, and for choosing vertices ut,vt on face 
ft we have 6 options. Observe that 


[SvTj 

X («")’■ < 

r=l 


L3vTj 

^(24)^" < (2d)i5^. 


As family Af, we output guarded separators X(5) for all the candidates S that are indeed a Voronoi 
separator; recall that this property can be checked in polynomial time, as well as the construction of 

Moreover, 
time and 


X(iS') takes polynomial time. Thus we have |AA| < (2d)^^^, which is exactly property ( 


Af can be enumerated by examining the candidates one by one within the requirec 


additional space. We are left with proving property (ii) 


Take any normal subfamily fF T V with IT"! = k, a nd le t M = Mj- and T-L = Aij be the 


corresponding Voronoi partition and diagram. By Lemma 4.5 we have that % is a. connected, 3 


regular multigraph (possibly with loops), \V{TL)\ = 2 | T'| — 4 = 2k —4, and \E{'H)\ = SlT "!—6 = 3k —6. 


4.17 


Since A: > 4, we have that \E{T-L)\ > 6 and Corollary 
noose 7 w.r.t. TL that has length at most y 4.5\V{TL)\ < 3y/k. By Lemma 


asserts that there exists a |-face-balanced 


4 . 1 f 


i) 


S ( 7 ) is a Voronoi 
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separator of the same length as 7 . Since C 2? is a normal subfamily, by Theorem 4.7 we have that 


4.21 


V{%) C £1, and hence 5 ( 7 ) uses only faces from £. This means that X(5(7)) G N. By Lemma 
we have that guarded separator 'K{S{j)) is |-interaction-balanced w.r.t. J^, and hence property (ii^ 
is proven. □ 


4.10 The algorithm 

We first briefly discuss the intuition. Suppose Z is a solution to the considered instance {G,'D,C,k), 
i.e., it is a normal family of exactly k objects from D. We can assume that k > 4, since otherwise 


the instance can be solved in polynomial time by brute force. Theorem 4.22 gives us a method 
to enumerate a small family of candidate guarded separators with a promise, that one of them 
separates the optimum solution evenly. More precisely, there is a guarded separator whose objects 
all belong to the optimum solution Z, and moreover after inferring all the information from this 
fact, i.e., excluding banned objects and covered clients, we arrive at the situation where every 
connected component of the remaining interaction graph contains only two thirds of the objects 
of Z. This gives rise to a natural Divide&Conquer algorithm: For every guarded separator X 


output by Theorem 4.22, we recurse into all the components of the modified interaction graph >C(X) 
for all parameters between 0 and and then combine the results using a standard knapsack 

dynamic programming. Note that in the new subinstances solved recursively we only modify the 
sets of objects and clients, whereas the underlying graph G remains unchanged. We now proceed to 
implementing this plan formally. 

Let us first introduce some notation. For an instance I = {G,T),C, k) of Disjoint Network 
Coverage, by Val[X] we denote the maximum revenue of a solution in X that should be reported by 
the algorithm. When X C X is a normal subfamily of objects, by Val[X, Z] we denote the revenue 
of Z in X. By N we denote the family of guarded separators computed for instance X using the 

Suppose X is some guarded separator in M. By cc(X) we denote the 


algorithm of Theorem 4.22 


set of connected components of X(X) = £ — (Q(X) U ban(X) U cov(X)), where £ is the interaction 
graph of instance X. For every G G cc(X), by X(C') and C{C) we denote the sets of objects and 
clients in C, respectively. Finally, we shall say that a vector of nonnegative integers k = (A:c)cecc(x) 
is compatible with 'K if kc < \k for each G G cc(X) and '^c&cc{X) kc = k — |Q(X)|. The set of 
vectors compatible with X will be denoted by comp(X). 

The following lemma is the main argument justifying the correctness of the algorithm. 

Lemma 4.23. Provided A: > 4, the following recursive formula holds for every instance {G, X, C, k) 
of Disjoint Network Coverage.- 


Val[G, X, C, fel = max max n(Q(X))-|- 

XeN kScomp(X) 


^ Ya.l[G,V{C),C{G),kc] 


(13) 


Cs ccC- 


Proof. Let X = {G,'D,C, k). For a connected component G on the right-hand side we shall denote 
Xc = {G,'D{G),C{G),kc)- Let L and R denote the left- and right-hand side of formula (13), 
respectively. 

We first show that L < R. If L = —00 then the inequality is obvious, so suppose that X C X 
is a normal subfamily of cardinality exactly k that maximizes the revenue on the left-hand side, 
i.e., Val[X] = Val[X, X]. By Theorem 4.22, there exists a guarded separator X G A/" that is 
|-interaction-balanced w.r.t. Z. In particular X is compatible with Z, so ban(X) n X = 0. 

For a connected component G G cc(X), let Zc = Z nX(C') and kc = \Zc\', observe that Zq is a 
normal subfamily of X(C'). Since ban(X) n X = 0, we infer that X is a disjoint union of Q(X) and 


sets Zc for G G cc(X). Hence X)cecc(x) = A: —|Q(X)|. Moreover, since ^ 


' is I-interaction-balanced 


47 












w.r.t. -H, we have that kc < for all C G cc(X). This implies that vector k := (^C')cecc(x) is 
compatible with X. 

Observe now that Val[X, Z] = n(Q(X)) + Val[Xc, Zc], where we denote Xc = 

{G,'D{C),C{C), kc)- Indeed, the cost of every object of Z is counted once on both sides of 
the equality, and the same holds also for the prizes of clients covered by Z: Prize of a client covered 
by an object of Q(X) is counted only in the term n(Q(X)) due to removing these clients when 
constructing >C(X), and every client not covered by Q(X) but covered by Z is counted only in the 
term Val[Xc,2^c] for the component C it belongs to. Recall also that Val[X] = Val[X, and 
Va^Xc] > Val[Xc',Xc] by the definition of Val[Xc]. Therefore, we have: 

L = Val[X] = Val[X,Z] =n(Q(X))+ Val[Xc, Xc] < n(Q(X)) + ^ Val[Xc] < R; 

Cecc(x) Cecc(x) 


the last inequality follows from the fact that pair (X, k) is considered in the maxima on the right-hand 
side of (13). 

Now we prove that L > R. For this, it suffices to show that 


Val[X] > n(Q(X)) + ^ Val[Xc] (14) 

Cecc(x) 


for every guarded separator X and every vector k = {kc)c&cc{X) such that X]cecc(x) kc = k — |Q(X)|; 
here, we denote Xc = {G,'D{C),C{C), kc)- Hence, let us fix such pair (X, k). For each C G cc(X), 
let Zc X T>{G) be a normal subfamily of cardinality kc that maximizes the revenue in the instance 
Xc; i.e., Val[Xc,Xc] = Val[Xc]. If such a family does not exist, then Val[Xc] = —oo and the 
claimed inequality is trivial. 

Let us now construct Z := Q(X) U Ucecc(x) ^c- Since all the sets in this union are pairwise 
disjoint, we have that \Z\ = k- 

We first verify that Z is normal. For the sake of contradiction, suppose that there are two objects 
G X that contradict the definition of normality. If p,p' G Q(X), then Q(X) would not be normal, 
a contradiction with the definition of a guarded separator. If p G Q(X) and p' ^ Q(X), then it would 
follow that p' G ban(X), which means that p' could not have been included into any Xc X T>(C), 
since these sets are disjoint with ban(X). If p,p' ^ Q(X) but p,p' GXc for the same component C, 
then we would have a contradiction with the normality of family Xc- Finally, if p G Xc and p' G Xc' 
for different components G,C' G cc(X), then we would have that pp' G X(X(X)), a contradiction 
with G and C' being different connected components of £(X). 

Now, we verify that 


Val[X,X] =n(Q(X))+ ^ Val[Xc,Xc]. (15) 

Cecc(x) 


Firstly, for the objects’ costs, since X = Q(X) Xc we have that the cost of every object of 

X is counted exactly once on each side of ([T§. Second, we check that the prize of every client (7 G C 
covered by X is counted exactly once in the right-hand side of ( |15[ ); note that this formula does 
not count the prize of any client not covered by X. If q is covered by Q(X), then 'K{q) is counted 
once in n(Q(X)) and in none of the terms Val[Xc,Xc], since q was removed when constructing the 
graph £(X). If now q is not covered by Q(X), then q belongs to some component G G cc(X). Since 
q is covered by X, there exists some p £ Z that covers g; observe that by the definition of £(X), 
for each such p we have that pq G E{C{X)), which implies that each such p must also belong to G- 
We infer that '^{q) is counted once in term Val[Xc,Xc], and is counted zero times in every term 
Val[Xc/,Xc/] for G' / G- This concludes the proof of formula (15). 
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Now, since Val[X] > Val[X,Z] and Val[Xc] 
obtain that 


Val[Xc,Xc] for every C G cc(X), by (15) we 


Val[X] > Val[X,X] =n(Q(X))+ Val[Xc, Xc] = n(Q(X)) + Val[Xc], 

Cecc(x) Cecc(x) 


which proves (14). Since X and k were chosen arbitrarily, we have thus proved that L > R. □ 


Lemma |4.23| justifies the correctness of the following recursive algorithm for computing the 
value of Val[G,X,C, k], summarized as Algorithm SolveDNC. As the border case, if fc < 3, then we 
compute the optimum revenue in a brute-force manner, by iterating through all the fe-tuples of 


the objects. Otherwise, we aim at computing Val[X] using formula (13). To this end, we run the 


algorithm of Theorem 4.22 to enumerate the family M of candidates for guarded separators separating 
evenly the solution. For each enumerated guarded separator X G AA we need to compute the value 
maxk6comp(x) n(Q(X)) + Y,ceccix)^^^iG,V{C),C{C),kc], since Lemma 


4.23 


computing Val[G, X, C, k] will boil down to taking the maximum of these values. 


asserts that then 
i'or a fixed guarded 

separator X, we investigate all the connected components of cc(X), and for each C G cc(X) and 
every integer £ with 0 < £ < |/c we compute the value Va^G, 'D{C),C{C),i] using a recursive call to 
Algorithm SolveDNC. Computing value maxkgcomp(x) n(Q(X)) -|- X]cecc(x) X(G), C(C'), kc] 

based on all the relevant values Val[G,X(G),C(G), £] amounts to running a standard knapsack 
dynamic programming algorithm that keeps track of the optimum revenue for selecting £ objects 
from the first i connected components of cc(X), for all £ G {0,1,..., A:} and i = 1, 2,..., |cc(X)|. 


Lemma 4.23 ensures that the algorithm correctly computes the value Val[X], so we are left with 
estimating the running time. 


Lemma 4.24. Algorithm SolveDNC, applied to an instance (G,'D,C,k) with |14(G)| = n, |X| = d, 
and \C\ = c, runs in time and space {dcn)^^^\ 

Proof. We shall say that a subcall (G,X',C', k') to Algorithm SolveDNC is trivial if V' = 0. Note 
that trivial subcalls are resolved by SolveDNC in polynomial time since family M is empty. 

Consider the application of SolveDNC to an instance {G,'D,C,k) with A: > 4, and let us count 
the number of subcalls to SolveDNC it invokes. We have that lA/"! < {2d)^^^, and for every X G AA 
we call SolveDNC on |cc(X)| • ([|A:J + 1) subinstances. At most d components of cc(X) contain an 
object of X, which means that at most d ■ ([|A:J + 1) < dk of these calls are non-trivial. In total, we 
have at most dk ■ { 2 d)^^^ < (2ci)^’^^ nontrivial subcalls. 

Let T* be the tree of nontrivial recursive subcalls of the application of SolveDNC to instance 
{G,'D,C, k). Observe that T* has depth 0(logk), since for each call of SolveDNC, in all the recursive 
subcalls the target size of the sought family is at most a |-fraction of the original one. For 0 < ^ < A;, 
let (f){£) be the maximum number of leaves in the subtrees of T* rooted at calls where the target 
size of the sought family is at most £ (or 1, if there are no such calls). Since in the recursive calls 
both the target family size and the total number of objects can only decrease, we have that (f>{£) 
satisfies that: 

(/){£) < (2d)^^^ • ^{\2£/3]) for 4 <£<k. (16) 

Moreover, since for £ < 3 Algorithm SolveDNC solves the instance by brute-force, we have that 


(/){£) = 1 for 0 < £ < 3. 


(17) 


Claim 4.25. For every function (/) : N — )• N that satisfies (16) and {If^, the following holds: 


<p{£) < (2(i)i7(3+vT)vX for0 <£<k. 


(18) 
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Input: An instance X = (G, P, C, /c) 
Output: Value Val[X] 


if A; < 3 then 

Compute the optimum revenue ret by iterating through all the /c-tuples of objects 

return ret 
end 

ret < -oo _ 

for X G M, enumerated using the algorithm of Theorem ^.22 do 
Compute n(Q(X)) 

Compute sets ban(X) and cov(X) 

Compute >C(X) = C — (Q(X) U ban(X) U cov(X)) 

Cl, C 2 ,... ,Cp e- Connected components of C{%) 
for i = 1 to p do 

for ^ = 0 to [|fcj do 


A[i][ll] ^ SolveDNC(G,2?(Ci),C(Cj),^) 


end 


end 


D[0][0] e- 0 

for £ = 1 to A: do 

I L»[0][£]^-oo 

end 

for i = 1 to p do 
for A* = 0 to A; do 

Z)[i][£] <-00 

for £' = 0 to min(£, [|A:J) do 

Z)[i][f] ^ma'x{D[i][£],A[i][i'] + D[i — l][£ — i']) 

end 

end 

end 

ret <— max(rei,n(Q(X)) + D[p][k — |Q(X)|]) 

end 

return ret 

Algorithm 1: Algorithm SolveDNC 


Proof. We prove the claim by induction. The base for £ < 3 is trivial by ( |17[ ), so let us show the 
inductive step for f > 4: 


m < (2d)i7-^.0(r2£/3i) 

= (2d)i7(i+\/6+2)vT ^ 


Here, in the first inequality we used (16) and in the second we used the induction hypothesis. 


Claim 


4.25 


implies that the subcall tree T* has at most < (2d)®^^ leaves. Since 

the depth of T* is C>(logA;), we infer that T* has C>(logA: • {2d)^'^'^) nodes. Now observe that for 
each call of SolveDNC, the amount of work used in the main body of this algorithm (excluding the 
recursive subcalls) is equal to the maximum possible number of guarded separators considered at 
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this step, i.e., times a factor polynomial in re, d, and c; here, we assume that the time 

used for resolving trivial subcalls is charged to the parent call. Hence, the total running time used 
by the algorithm is at most • (cre)^^^\ as requested. To argue that 

the algorithm runs in polynomial space, observe that at each moment the algorithm keeps track of 


the data stored for 0(logA:) recursive calls of SolveDNC, and, due to the algorithm of Theorem 4.22 
using only polynomial working space, each call uses only polynomial space for its internal data. □ 


Thus, Lemmas 4.23|and 4.24 conclude the proof of Theorem |3.1 


5 Hardness results 


In this section, we prove the lower bounds in Theorems 1.8| - 1.11 suggesting that for many natural 
covering problems the time brute force algorithms are almost optimal. Let us start with a 

very simple reduction showing such a lower bound in the setting of planar graphs. Patra§cu and 
Williams gave a very tight lower bound for Dominating Set, assuming the Strong Exponential 
Time Hypothesis (SETH). 


Theorem 5.1 (|44j). Assuming SETH, there is no f{k) ■ P time algorithm for Dominating 
Set, where re is the number of vertices of the input graph G. 


Using this result, the proof of Theorem 1.8 is very simple and transparent. 


Restated Theorem 1.8 (covering vertices with connected sets, lower bound). Let G be 

a planar graph and let V be a set of connected vertex sets of G. Assuming SETH, there is no 
f{k) • {\V\ + |U(G)|)^“^ time algorithm for any computable function f and any e > 0 that decides if 
there are k sets in V whose union covers |U(G)|. 


Proof. Let H be an arbitrary graph with re vertices ui, ..., Vn- Let G be a star with leaves ui, 
..., Vn and center x. For every vertex Vi G V{H), let us introduce a set Si into V that contains 
Vi, every neighbor of Vi in H, and the center vertex x of G; clearly, the set Si is connected in G. 
It is easy to see that k vertices of H form a dominating set of G if and only if the corresponding 
k sets of H cover every vertex of G. Therefore, an f{k) • (|X1| + |U(G)|)^“'^ time algorithm for the 
covering problem on G would give an f{k)n^~^ time algorithm for solving Dominating Set on H] 
this contradicts SETH by Theorem |5.1[ □ 


[mi 


In Sections 5.1-5.31 we give similar lower bounds for geometric objects, proving Theorems 1.9 


5.1 Convex polygons 

The proof of Theorem |1.9| is again a simple reduction from Dominating Set. We observe that 
if a family of points lies on a convex curve, then any subset of them is in convex position, hence 
convex polygons can cover arbitrary subsets of these points. This simple idea was also used, for 
example, by Har-Peled 123, but here we use it to obtain a different conclusion (tight lower bound 
on the exponent instead of APX-hardness). 

Restated Theorem 1.9 (covering points with convex polygons, lower bound). Let H be 

a set of convex polygons and let V be a set of points in the plane. Assuming SETH, there is no 
f{k) ■ {\V\ + 17^1)^“*^ time algorithm for any computable function f and e > 0 that decides if there 
are k polygons in H that together cover V. 
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P8 = ((n + l)2,(n+l)2) 


Po 



= W 


Figure 13: Proof of Theorem 1.9 The points on a parabola and the convex polygon corresponding 
to a vertex whose closed neighborhood is {^ 2 ,^ 5 }. 


Proof. Given an instance {G,k) of Dominating Set with n = |F(G)|, we construct a set P of n 
convex polygons and a set V of n + 2 points such that G has a dominating set of size k if and only 
if V can be fully covered by k convex polygons from P. Therefore, an f{k) ■ n^~'" time algorithm 
for the problem would imply an algorithm for DOMINATING Set with similar running time, which 


would contradict SETH by Theorem 5.1 


We define the point set V = {pi = (i(n + 1), | 0 < i < n + 1}, which is a set of n + 1 points 

on a parabolic curve (see Figure [T^. Let ui, ..., be the vertices of G. For every 1 < i < n, we 
introduce into P a convex polygon corresponding to Uj. Suppose that the closed neighborhood of Vi 
is Vj^, ..., Vj^ for some 1 < ji < j 2 < • • • < id < n. Then Vi is represented by a convex polygon 
Pi with vertices (0, (n + 1)^), po, pj^, Pj^, • • ■, Pj^, Pn+i- As the parabola is a convex curve, this 
polygon Pi is convex and it covers exactly the points po, pj.^, pj.^^ • • •) Pjr^ Pn+i- That is, if vj is not 
in the closed neighborhood of Vi, then pj is outside Pj. It is easy to see that a subset of P covers V 
if and only if the union of the closed neighborhoods of the corresponding vertices of G cover V{G), 
or in other words, the corresponding vertices of G form a dominating set. □ 


5.2 Thin rectangles 

In the proof of Theorem jl.lOt we are reducing from the Partitioned Biclique problem: given 
a graph G with a partition {A\,..., A^, Bi,..., B/^) of the vertices, the task is to find 2k vertices 
ai G Ai, ..., Ok G Ak, bi G Bi, ..., bk G Bk such that a* and bj are adjacent for every 1 < i, j < k. 
There is a very simple reduction from Clique to Partitioned Biclique (see, e.g., m) where the 
output parameter equals the input one, hence the lower bound of Chen et al. [8] can be transferred 
to Partitioned Biclique. 

Theorem 5.2. Assuming ETH, Partitioned Biclique cannot be solved in time f{k)nA^'> for 
any computable function f. 

The proof of Theorem |1.10| is a parameterized reduction from Partitioned Biclique. 

Restated Theorem 1.10 (covering points with thin rectangles, lower bound). Consider 
the problem of covering a set V of points by selecting k axis-parallel rectangles from a set P. 
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Assuming ETH, there is no algorithm for this problem with running time f{k) ■ {\V\ + for 

any computable function f, even if each rectangle in V is of size 1 x k or k x 1. 


Proof. We prove the theorem by a reduction from Partitioned Biclique. Let G be a graph with 
a partition (^i,..., Ak, Bi,..., Bk) of the vertices. Without loss of generality, we assume that every 
class of the partition has the same size n. For 1 < i < A:, we denote the vertices of Ai and Bi by aij 
and bij (1 < j < n), respectively. 

Construction. We construct a set TZ of axis-parallel rectangles and a set V of points such that 
there are k' := 8k rectangles in TZ covering 'P if and only if G contains the required partitioned 
biclique. For simplicity of notation, we construct open rectangles, that is, selecting a rectangle does 
not cover the points on its boundary. However, it is easy to modify the reduction (by decreasing the 
size of each rectangle slightly) so that it works for closed rectangles as well. 

First, for every 1 < i < A, we define the rectangles (see Figure [Td]) 

= (4i-0.5,4i-b0.5) x {4:{k + 1), 8{k + 1)) 

= (4i - 0.5,4i -h 0.5) x {-A{k -h 1), 0) 

= (4i-l,4i) X (0,4(A: + 1)) 

Vf- = {U,Ai + l) X (0,4(A + 1)) 

and for every 1 < j < A, we define the rectangles 

Hj = (0,4(A: + 1)) X (4j,4j + l) 

Ff = (0,4(A + l))x(4j-l,4j) 

= {-4{k + 1), 0) X (4j - 0.5, Aj + 0.5) 

Hf = (4(A; -h 1), 8{k -b 1)) x (4j - 0.5,4j -b 0.5) 

These rectangles themselves do not appear in the set TZ, but they will used be for the definition 
of rectangles in TZ, as follows. For a rectangle R, we use the notation R + (x, y) to denote the 
rectangle obtained from R by shifting it horizontally by x and vertically by y. We use the notation 
\ etc. for the diagonal vectors (1,1), (1, —1) etc. with unit projections; then, e.g., A- \ denotes 
the vector (A, — A). 

Let e = l/(100n). For every 1 <i <k and 1 < a < n, we introduce the following rectangles 
into TZ (see the directions shown by the arrows in Figure [l^: 

+ ae- \ 

v.l = yl' + \ 

Furthermore, for every 1 < j <k and 1 < /I < n, we introduce the following rectangles: 

= Hf + 0,- \ 

^ 

This completes the description of TZ] note that \TZ\ = 8kn. To define "P, for every 1 < f < A: and 
1 < a < n, we introduce the points (see again the directions shown by arrows in Figure [T^ 

vfjf = (4i - 0.5, 4{k -b 1)) + (a + 0.5)e- Z' u™ = (4i -b 0.5, 4(A: -b 1)) + {a + 0.5)e- \ 

= (4^ - 0.5, 0) + (a + 0.5)e- \ = (4i + 0.5,4(A: + 1)) + (a + 0.5)e- ^ 
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Figure 14: Proof of Theorem 


1.10 


The rectangles on the left and the right are shortened to save 


space. The directions of the arrows on the rectangles show the offsets of the rectangles etc. as 
a and /3 grow. Similarly, the arrows at etc. show the offset of the vertices hj^. 
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and for every 1 < j < k and 1 < /3 < n, we introduce the points 

hl^ = (0,4j + 0.5) + (/? + 0.5)e- / hJJ = {4{k + 1), 4j + 0.5) + (/? + 0.5e)- \ 

hfj = (0,4j - 0.5) + {/3 + 0.5)e- \ hfj = {4{k + 1), 4j - 0.5) + (/3 + 0.5)e- ^ 

We will use the notation ■= {vf^ : 1 < a < n} etc. for these sets of points. Notice that these 

points are not on the boundary of any of the rectangles (because of the terms 0.5e), thus it does not 
matter for the covering of these points that the rectangles are open. 

Finally, for every 1 < i, j < k, we add a set Pi^ of points to V the following way: for every 

1 < a, /3 < n, the point (4z + ae, 4j + f3e) is in Pij if ai^a and are not adjacent in G. This 

completes the description of construction, we have 


r = [J{v, 


ii)L yTFt ^BFt\ 


i=l 


uU(^' 


h,TL y ^h,TR y ^h,BL y 


u 


u 


p 




1<2 J</C 


Note that \V\ < 8kn + /c^n^. 

We claim that we can cover the points in V with 8k rectangles of TZ if and only if the Partitioned 
Biclique instance has a solution. The intuition behind the construction is the following. A solution 
of size 8k has to select exactly one rectangle of each of the 8k types ®tc. The fact that 


the point set U U p^P^ jjv,br covered by the four selected rectangles of 

the form ensures that the selection of these rectangles are synchronized, 

that is, we have to select some 1 < a* < n. Similarly, for every I < j < k, 

we have to select some 1 < j3j < n. The role of the point set Pij is 

to ensure that vertices Oi^ai and are adjacent in G. Indeed, if they are not adjacent, then the 
point (4i + ae, 4j + /3e) is in Pij, but it is not covered by any rectangle: it is on the boundary of 
each ofP4,P4,//J^^.,P4.’ 

For the formal proof of the equivalence, let us start with a useful preliminary observation. 

Claim 5.3. For every I < j < k and 1 < aT,aB,oiL,aR < n, the rectangles 

cover all the points in U U p^P^ y p^PR = q,^. 


CtL’ 


Proof. Assume first that ax = oib = oil = oir = a*. Consider a point vj^. If a < a*, then 
rectangle covers this point: the vertical coordinate of this point is 4{k + 1) + ae + e/2 and 
the top boundary of is at 4(A: + 1) + a*e. If a > a*, then rectangle P^^* covers the point: 
the horizontal coordinate of this point is 4i — 0.5 + (a + 0.5)e and the left boundary of is at 
4i — 0.5 + a*e. In a similar way, we can verify that the four rectangles cover vf^, vfF: vfF fo^' 
every 1 < a < n. 

Suppose now that the four rectangles cover all the 4n points. If a^ < cut, then point is 
not covered: rectangle Pj^^ covers only the points vj^, ..., of while rectangle 

covers only the points ..., Uyf of Therefore, we have ai > cxt- In a similar way, we 

can prove ax > o:r > aB > cil P ctx, implying that all these values are equal. j 


We can prove an analogous statement for the vertical rectangles. 

Claim 5.4. For every 1 < j < k and 1 < fdx, (^b, h, I3 r < n, the rectangles 

Pfpg cover all the points in p^P^ y p^P^ y ph,BL ^ ph,BR if jdx = Pb = Pl = Pr- 

Biclique covering rectangles. Suppose that vertices ai^ai, ■ ■ ■, Ofc.ofe, • • •) form 
a biclique in G. We claim that the 8k rectangles 
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V-^ V- 

t,ai > t,a. 


^ for 1 < z < /c, and 

forl<J<A: 


cover every point in V. By Claim 


5.3 


for every 1 < z < A:, the rectangles 


5.4 


1,0-1^ Z^OCi^ " Z,OCi^ " 

for every 1 < J < fe, the 


cover every point in u U u p^^^R Claim 

rectangles HJ^., every point in u k jyh,BL ph,BR^ Consider 

now a point p = (4z + ae, 4j + /3e) G Pij. If a < Ofj, then V/"^. covers p; if a > a*, then covers p. 
If /3 < 13j, then covers p; if /3 > /3j, then I/?g, covers p. Therefore, we have a problem only if 
a = ai and {3 = (3j. However, we know that vertices and bj^p. are adjacent in G, which means 
that the point p = (4z + cte, 4j + /3e) was not added to Pij. Therefore, the selected 8k rectangles 
indeed cover every point in P. 

Covering rectangles =► biclique. For the proof of the reverse direction, suppose that V was 
covered by 8k rectangles from P. For every 1 < z < /c, at least one rectangle of the from has 
to be selected: no rectangle of some other form can cover the point vj^ with vertical coordinate 
4(A: + 1) + ne + e/2. Similarly, at least one rectangle of each of the forms 4/^, 4/^ has to be 

selected. An analogous statement holds for every 1 < j < A: and rectangles of the form 

Hfj^- As exactly 8k rectangles were selected, exactly one rectangle has to be selected for each 


type. Consider some 1 < z < A: and suppose that , 4/^^, 4/^^, Hj were selected. As these 


rectangles have to cover Claim 5.3 implies that ar = cxb = ai = otR', 






were selected. 

J^Pj 


let us define a* to be this number. Similarly, using Claim [5141 we can show that for every 1 < j < A:, 
there is a 1 < /!,• < n such that the rectangles Hja, 

J'lPj . „ „ 

We claim that ai^ai, ■ ■ ■, <^k,a^i ^i,Pii ■ ■ ■ ■, form a solution to the Partitioned Biclique 
problem. Suppose for a contradiction that ai^ai and bj^a- are not adjacent in G. Then the point 
(4z + ajC, 4j + (3je) is in Pij. We arrive to a contradiction by showing that this point is not covered 
by any of the selected rectangles. Indeed, this point is on the boundary of each of 
Hfa, and recall that all the rectangles are open. Therefore, we have shown that the Partitioned 
Biclique instance has a solution. □ 


5.3 Almost squares 


The goal of this section is to prove Theorem l.Il which is another lower bound for covering points 
with axis-parallel rectangles. In Theorem 1.10, we had only two types of “thin” rectangles. In this 
section, we consider sets of rectangles of many different sizes, but they are all “almost unit squares”: 
every rectangle has width and height in the range [1 — eo, 1 + eo] for some fixed eo > 0. 

Restricting the problem to almost unit squares makes the hardness proof significantly more 
challenging. Let us observe that the proof of Theorem 1.10 crucially relied on the fact that the 


selection of the horizontal rectangles interacted with the selection of the vertical rectangles 

Pj^ for every 1 < z, j < A:. Intuitively speaking, the selection of these 0{k) rectangles were 
constrained by 0{P) interactions. The reason why such a large number of interactions could be 
reached is because a wide horizontal rectangle could be intersected by 0{k) independent thin vertical 
rectangles. However, if every rectangle is almost a unit square, then this is no longer possible. After 
trying some possible conhgurations that may be useful for a hardness proof, one gets the impression 
that an almost unit square can have independent interactions only with 0(1) other almost unit 
squares. This means that we can have a total of 0{k) interactions, which does not seem to be 
sufficient to express a problem such as Clique or Partitioned Biclique, as these problems have 
a richer interaction structure with 0{P) interactions. Therefore, it seems difficult to reduce from 
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these problems with only a linear blowup of the parameter, which is necessary for the lower bounds 
that we want to prove. 

To get around these limitations, we are reducing from a problem where the interaction structure 
is sparser and a lower bound ruling out (or so) algorithms is known even with only 0{k) 

constraints. The Partitioned Subgraph Isomorphism problem is similar to Partitioned 
Biclique, but now the graph we are looking for can be arbitrary, it is not necessarily a biclique. 
The input consists of a graph H with vertex set {ui ,... ,Uk} and a graph G whose vertex set is 
partitioned into k classes Vi, ..., I 4 . The task is to find a mapping /r : V{H) —>• V{G) such that 
/i(ui) G Vi for every 1 <i <k and is a subgraph embedding, that is, if Ui and Uj are adjacent in H, 
then /u(ui) and nju j) ar e adjacent in G. This problem contains Partitioned Biclique as a special 

implies that it cannot be solved in time f{k)n°^^^ for any computable function 


5.2 


case, thus Theorem 

/. However, if we express Partitioned Biclique with Partitioned Subgraph Isomorphism, 
then H has 0{k‘^) edges. Interestingly, the lower bound remains valid (up to a log factor) even if we 
restrict H to be sparse. 

Theorem 5.5 ()39j). Assuming ETH, Partitioned Subgraph Isomorphism cannot he solved 
in time (where k = \ V{H)\) for any computable function f, even when H is restricted 

to be a 3-regular bipartite graph. 

We prove a lower bound first for an intermediary problem, covering points by 2-track intervals, 
and then obtain Theorem l.Il using a simple reduction of Chan and Grant [7]. Many natural 


optimization problems on intervals (such as finding a maximum set of disjoint intervals or covering 
points on the line by the minimum number of intervals) are polynomial-time solvable. One can 
obtain generalizations of these problems by considering c-intervals instead of intervals, that is, the 
objects given in the input are unions of c intervals and we have to select k such objects subject to 
some disjointness or covering condition [lEniEIlEaEIllsS]. These problems are typically much 
harder than the analogous problems on intervals. A special case of a c-interval is a c-track interval: 
we imagine the problem to be defined on c independent lines (tracks) and each c-track interval 
consists of an interval on each of the tracks. Many of the hardness results for c-intervals remain 
valid also for c-track intervals. 

Let us formally define the intermediary problem and prove the lower bound we need. A 2-track 
interval is an ordered pair (J^,/^) of intervals on the real line; for concreteness the reader may 
assume that the intervals are closed, but this does not have any influence on the validity of the 
claims to follow. We say that 2-track interval covers a point x on the first (resp., second) 

track if X £ (resp., x G /^). A 2-track point set V is an ordered pair (P^jP^) of two sets of real 
numbers. We say that that a set X of 2-track intervals covers a 2-track point set (P^,P^) if 

• for every x G P^, there is a 2-track interval (I^,/^) G X that covers x on the first track (i.e., 
X G I^) and 

• for every x G P^, there is a 2-track interval (I^, /^) G X that covers x on the second track (i.e., 
X G P). 


We prove a lower bound for the problem of covering a 2-track point set by selecting k 2-track intervals. 
Note that this problem is related to, but not the same as, the Dominating Set problem for 2-track 
interval graphs, for which W[l]-hardness results are known |2n( 132] (however no ETH-based lower 
bound was stated). 

Theorem 5.6. Consider the problem of covering a 2-track point set V = (P^,P^) by selecting k 
2-track intervals from a set X. Assuming ETH, there is no algorithm for this problem with running 
time f{k) ■ (|P| + \X\YVn^&k) any computable function f. 
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[Ani + j, Ani + n + j] [Ani + 2n + j, 3n + j] 

[Ani, Ani + j] ^ [Ani + n + j, Ani + 2n + j] [Ani + 3n + j, 4n] 



Ani Ani + n Ani + 2n Ani + 3n Ani + 4n 


Figure 15: Proof of Theorem 


5.6 


Let n = 8 and j = 3. The figure shows the 2-track intervals 
. Notice that together they cover [4ni, 4ni -|- 4n] on the second 


(light gray) and Iij ^2 (dark gray) 
track. On the hrst track, the three intervals of length n (shown by red, green, and blue) are needed 
to complete the covering of the range [Ani, Ani + An\. These intervals can be provided by the 2-track 


intervals F 


. 0 ) 


el ? 


1*2 1 where e'^ is an edge of G incident to Vi^j and having color c. 


Proof. We prove the theorem by a reduction from Partitioned Subgraph Isomorphism. Let 
H and G be two graphs, let V{H) = {ui,..., Uk}, and let {Vi ,..., 14) be a partition of V(G). By 
copying vertices if necessary, we may assume that every V) has the same size n; let us denote by 
..., Vi^n} the vertices in V). As Lf is bipartite and 3-regular, the edges of H are 3-colorable. Let 
us fix a 3-coloring of the edges of H with colors in {0,1, 2}. If there is an edge Ui^Ui^ in H having 
color c G {0,1, 2}, then we will also refer to the edges of G between and as having color c. 
We may assume that every edge of G receives a color this way, that is, there is no edge induced by 
any V) and there can be an edge between and Vi^ only if Ui^Ui^ is an edge of H (every other 
edge of G is useless for the purposes of this problem). As H is 3-regular, the two partite classes of 
H have to be of the same size; we may assume without loss of generality that {ui,..., Uh/ 2 } 
{ttfc/ 2 +i,..., Uk} are these two classes. 

Construction. We define the set I of 2-track intervals the following way (see Figure [T^. First, 
for every 1 < i < k/2 and 1 < j < n, let us introduce the following two 2 -track intervals into I: 

• h,j,i '■= {[Ani, Ani + j], [Ani, Ani + n — j]) 

• '■= {[Ani -|- 3n -|- j, Ani + An], [Ani + n — j, Ani + 4n]) 

For every k/2 + 1 < i < k and 1 < j < n, we proceed similarly, but exchanging the role of the two 
tracks. That is, we introduce the following two 2-track intervals into I: 

• ■= {[Ani, Ani + n — j], [Ani, Ani -|- j]) 

• Ii,j ,2 '■= {[Ani + n — j, Ani + An], [Ani -|- 3n -|- j, Ani + 4n]). 

Next, we introduce a 2-track interval into I for every edge of G. Let e = be an edge of G 

with 1 < ia < fc/2 and k/2 + 1 < it < k. Suppose that this edge e has color c G {0,1, 2} in the edge 
coloring we fixed. Then we introduce 

• Ig := {[Ania + cn + ja, Ania (c -h l)n -f ja], [Anib + cn + jt, Anib -k (c l)n -k jb]) 

into the set I. This completes the construction of I; note that \I[ = 3.5nk (as H is 3-regular). 

Let V = {4n -k 0.5,4n -k 1.5,..., 4n(A: + 1) — 0.5}. Let k' = 2[V{H)[ -k [E{H)[ = 3.5fc. We claim 
that it is possible to select k' members of X to cover the 2 -track point set {V,V) if and only if the 
Partitioned Subgraph Isomorphism instance has a solution. 
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Subgraph embedding 2-track intervals. Suppose first that vertices ..., form 
a solution of the Partitioned Subgraph Isomorphism instance. Then we can cover {V,V) by 
selecting the following set X' of 2-track intervals: 


• and Ii^si ,2 for 1 < i < k, and 

• for every edge of H (with 1 < ia < k/2 and k/2 + 1 < it < k), edge e = Vi^^si^Vi^^si 

has to be an edge of G; we select the corresponding 2 -track interval /*. 


We claim that every point in [4n, 4n(A; -|- 1)] is covered by X' on both tracks; this clearly implies that 
{V,V) is covered. Consider the range [4m, 4ni-|-4n] for some 1 < f < k/2 (the case A:/2-|-l < i < k is 
similar). On the second track, covers [4m,4m-|-n —Sj] and covers [4m-|-n — Sj, 4m-|-4n], 
thus the range is indeed covered (see Figure 15). Let k/2 < < k he the three neighbors 

of Ui in H, with the edge UiUic having color c. This means that the edge = Vi^SiVi‘^,Sic exists in 
G and I/c was selected into X' for every c G {0,1,2}. Recall from the definition that I/c covers 
[Ani + cn + Sj, 4ni -|- (c -|- l)n -|- Sj] on the first track. Now we have that, on the first track, 


• [4nf, 4m -|- Sj] is covered by 

• [Ani + Si , Ani -)- n + Sj] is covered by /*o, 

• [Ani + n + Si, Ani -|- 2n -|- Sj] is covered by I*i, 

• [Ani -|- 2n -|- Sj, 4ni -|- 3n -|- Sj] is covered by /* 2 , and 

• [Ani -|- 3n -|- Sj, 4ni -|- 4n] is covered by Ii,j, 2 - 

Thus [4m, Ani + An] is indeed covered and we have shown that {V,V) is covered by X'. 

2-track intervals subgraph embedding. Suppose now that there is subset X' XX of size 
k' that covers the 2-track point set {V,V). Let us define the following subsets of X: 


• := {lij^i I 1 < j < nj, 

• := {Ii,j,2 I 1 < J < nj, 

• X/^ ■= {1/ I edge e has color c and has an endpoint in Vi}. 

Notice that if Ui^Ui^ is an edge of H having color c, then X*^ ^ = X*^ ^ (as every edge with color c 
and having an endpoint in has its other endpoint in V)^). Looking at the definitions, we can 
observe that 


• for 1 < i < k, point Ani + 0.5 on the first track is covered only by members of 

• for 1 < i < k/2, point 4ni -|- 4re — 0.5 on the second track is covered only by members of Xi^ 2 , 

• for k/2 + 1 < i < k, point 4ni -|- 4n — 0.5 on the first track is covered only by members of Xj_ 2 , 

• for 1 < f < k/2, 0 < c < 2, point Ani -|- (c -|- l)n -|- 0.5 on the first track is covered only by 
members of I*^, and 

• for A;/2-|-l<f<fc, 0<c<2, point Ani -|- (c -|- l)n -|- 0.5 on the second track is covered only 
by members of I*^, and 

That is, at least one member has to be selected from each of the 5k sets Xj^i, Xj_ 2 ; X*^. Note that 
every member of X is contained either in exactly one of these sets (in Xi^i or Xi ^2 for some 1 < i < k) 
or in exactly two of these sets (in X*^ ^ and X*^ ^ for some 1 < ia < k/2, k/2 + 1 < it < k, c £ { 0 , 1 , 2 }). 
Therefore, the only way to cover the 5k points enumerated above by selecting k' = 3.5A: members of 
X is to select exactly one member of each Xj^i, Xi^ 2 , and X*^. 

Suppose that for some f < i < k/2, the solution X' contains the unique members lij £ Xjq, 
Iij> £ Xj^2, and /*c G X*^ for c G {0,1,2}. By definition of X*^, edge has an endpoint Vijc £ Vi. 
We claim that j = j' = = j^. Recall that covers [Ani, Ani + n — j] on the second 

track, while /ij ',2 covers [Ani + n — j',Ani + 4n] on the second track. Thus if j > j' , then point 
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Ani + n — j' — 0.5 would not be covered on the second track; hence we have j' > j. Observe next 
that covers [4m,4ni + j] and /*o covers [Ani + , Ani + n + j^] on the first track. Thus we infer 

that j > otherwise, point Ani + — 0.5 on the first track would not be covered. In a similar 

way, we can show that and j'^ > j' also hold. Therefore, we get the chain of inequalities 

j > f > / > f > j' > j, 


implying that we have equalities throughout, as claimed. 

In the previous paragraph, we have shown that for every 1 < i < k/2, there is a 1 < Sj < n 
such that X' contains the unique members li^Si £ Ii,Si £ Ii, 2 , and I*c £ T*^ for c G {0,1,2} 
with Vi^si being an endpoint of e^. In a similar way, we can prove the same statement for every 
k/2 + 1 < i < k. We claim that ..., Vk^sj^ form a solution. Suppose that for some 

A < ia < k/2, k/2 + 1 < < A: is an edge of H, having color c G {0,1,2}. We have to show 
that and are adjacent in G. Recall that I*^ ^ ^ and consider the unique member 

/* G X* ^ selected into X'. As shown above, this means that edge e of G has an endpoint £ Via 

and also an endpoint G what we had to show. □ 


We are now able to prove Theorem 1.11 by a simple and transparent reduction from covering 
points by 2-track intervals. The basic idea is that if we have two parallel lines, then for any pair 
of segments on these two lines, there is a rectangle that intersects these lines exactly at these two 
segments. Therefore, by arranging the points of a 2-track point set on two parallel lines, we can 
easily turn the problem of covering by 2 -track intervals into a problem of covering by rectangles. 
Chan and Grant [7] used the same reduction to establish the APX-hardness of the problem (by 


reducing a certain special case of covering by 2-track intervals). Thus once we have Theorem 5.6 


our lower bound for covering by 2 -track intervals, the reduction of Chan and Grant [7j implies 
Theorem II.Ill 


Restated Theorem 1.11 (covering points with almost squares, lower bound). Consider 
the problem of covering a set V of points by selecting k axis-parallel rectangles from a set V. 
Assuming ETH, for every eo > 0, there is no algorithm for this problem with running time 
f{k) ■ {\V\ -|- for any computable function f, even if each rectangle in V has both 

width and height in the range [1 — eo, 1 + eo]. 


Proof. The proof is by reduction from the problem of covering points by 2-track intervals. Let X 
be a set of 2-track intervals, let V = ,V‘^) be a 2-track point set, and let k be an integer. By 

rescaling if necessary, we may assume without loss of generality that every point in and is in 
the interval [0,1] and hence every (J^,/^) G X satisfies 1^,1^ C [0, Ij. As shown in Figure 


16 


we 


arrange the points of (resp., R^) on the line y = x + 1 (resp., y = x — 1). That is, we construct 
a 2-dimensional point set V* that contains 


• point (eos, 1 + eos) for every s £ and 

• point (1 -t- cqs, eos) for every s £ . 


Then we represent every 2-track interval ([ai, 6 i], [02, 62]) G X by a rectangle whose horizontal 
projection is [eoui, 1 -|- 6062] and vertical projection is [6902 > 1 + eo^i]; note that in this proof we 
assume that the rectangles are closed, but to obtain the same result for open rectangles we can 
simply make them slightly larger. Observe that the sizes of both projections are between 1 — eo and 
1 -|- cq. Moreover, these rectangles faithfully represent the way the 2-track intervals cover the points 
in For example, 2-track interval ([ai, 6 i], [ 02 , 62 ]) covers point s on the hrst track if and 

only if oi < s < 5i and the corresponding rectangle covers point (cqs, 1 + cqs) exactly under the 
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1 + eo^i 



1 + £002 



+ 


Figure 16: Proof of Theorem 


1.11 


The rectangle [eooi, 1 + £ 0 ^ 2 ] x [£ 0 ^ 2 ,1 + £ 0 ^ 1 ] covers exactly 


those points that correspond to points in [m, 61 ] on the first track and or to points in [ 02 , ^> 2 ] on the 
second track. 


same condition (if s < ai, then the point is to the left of the rectangle; if s > 61 , then the point is 
above the rectangle). Therefore, we can cover the point set V* with k of the constructed rectangles 
if and only if we can cover the 2-track point set V with k of the 2-track intervals from I. □ 
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